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The purpose of the review - to provide affordable to a wider audience a description of kinematics 
of the cosmological expansion, and dynamic interpretation of this process. The focus will be the 
accelerated expansion of the Universe. Will be considered by virtually all major opportunity to 
explain the accelerated expansion of the Universe, including an introduction to the energetic budget 
of the universe of so-called dark energy, a modification of Einstein's equations, and a new approach 
based on the holographic principle. 



I. INTRODUCTION 
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There are two fundamental reasons making the question put in the title so principal. Firstly, in common sense, the 
correct answer is necessary (though insufficient) to predict the ultimate fate of the Universe. To be more specific, one 
must know current values of the kinematic parameters (velocity, acceleration etc.) in order to fix initial conditions 
necessary for solution of the differential equations describing dynamics of the Universe. Secondly, description of either 
decelerated or accelerated expansion of the Universe requires absolutely distinct cosmological models. Conventional 
substances (non-relativistic matter and radiation) from one hand and general relativity from the other explain the 
(-H ■ decelerated expansion. In order to obtain the accelerated expansion one must either change the composition of the 
^ ^, Universe rather radically, or make even more responsible conclusion: the fundamental physical theories lying in the 
basis of our understanding of the World, are essentially wrong. Common experience tells us that it is very unlikely 
to answer the posed question in near future. Already 80 years passed after the discovery of the Universe expansion, 
and we still cannot establish the source responsible for the initial velocity distribution. "The Big Bang" is nothing 
^ \ but euphemism used to hide limitation of our knowledge. What does actually expand? This question is still sharply 
controversial: "How can vacuum expand?" Q. Even more miraculous, how can vacuum expand with acceleration 
[ or deceleration? "...how is it possible for space, which is utterly empty, to expand? How can nothing expand? The 
, answer is: space does not expand. Cosmologists sometimes talk about expanding space, but they should know better." 

, No generation of homo sapiens can vanquish the temptation to believe that it is chosen to acquire complete and 
' final understanding of Nature and Universe. In particular, the first half of the last century, with discovery of Universe 
\ expansion by Hubble and due to now fundamental physical theories, relativity and quantum mechanics first of all, was 
QQ ' marked by new cosmological paradigm — the Big Bang model. The main basis of the model is the celebrated Hubble 
law, discovered in twenties of the last century and supported by all types of up-to-date cosmological observations. 
. The Big Bang model was able to give satisfactory explanation for thermal evolution of cooling Universe, predicted 
' existence of relict radiation, correctly described relative abundance of light elements and many other features of the 
^ , Universe. 

\^ ' Towards the end of the last century it was commonly hoped that the Big Bang model supplemented by the Inflation 
rS [ Theory represented an adequate (at least as a first approximation) model of Universe. However this hope was never 
■ to come true. The cosmological paradigm had to change for the reason of observations performed with ever increasing 
precision. Right after the Hubble discovery the cosmologists tempted to measure the deceleration of the expansion, 
attributed to gravitation. They were so firmly confident to discover exactly above stated effect, that the corresponding 
observable was named the deceleration parameter. However, in the year 1998 two independent collaborations 0,13], 
studying distant supernovae, presented convincing evidence for the fact that the Universe expansion is accelerated. It 
turned out that the brightness decreases in average much faster than it was commonly believed according to the Big 
Bang model. Such additional dimming means that the distance, attributed to given redshift, was somewhat underes- 
timated. And it in turn means that the universe expansion accelerates: in the past the Universe expanded slower than 
nowadays. The discovery of cosmological acceleration is probably one of the most important observations not only in 
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modern cosmology but in physics in general. Accelerated expansion of Universe the most directly demonstrates the 
fact that our fundamental theories are either incomplete, or, even worse, misleading ]. 

Physical origin of the cosmological acceleration still remains greatest miracle. As we already mentioned above, if 
the Universe is solely composed of matter and radiation, it must decelerate the expansion. If the expansion in fact 
accelerates, we have to choices, each of which forces us to revise our basic physical concepts. 

1. Up to 75% of energy density in the Universe exists in form of unknown substance (commonly called the dark 
energy) with high negative pressure, providing the accelerated expansion. 

2. General Relativity theory must be revised on cosmological scales. 

We should remark that, besides the two above cited radical possibilities to put the theory in agreement with the 
observations, there still exist an obvious conservative way to get rid of the problem: it implies more adequate utilization 
of available theoretical possibilities. D.Wiltshire expressed this idea in the following way: I will take the viewpoint 
that rather than adding further epicycles to the gravitational action ^ the cosmological observations which we currently 
interpret in terms of dark energy are inviting us to think more deeply about the foundations of general relativity. The 
above mentioned term " more adequate utilization of the available possibilities" need an explanation. Let us give an 
example. In the context of cosmological models based on homogeneous and isotropic Universe, in order to explain the 
observed acceleration one have to involve a new form of matter with negative pressure — the dark energy. According 
to an alternative interpretation [l^ - fTBj . the accelerated expansion of the Universe follows from deviations from the 
homogeneity. It is assumed that mass density in the Universe is considerably inhomogeneous on the scales smaller 
than the Hubble radius. In order to transit to effective description of homogeneous and isotropic Universe, one should 
average and/or smooth out the inhomogeneities up to some properly chosen scale of the averaging. In such an averaged 
Universe one can define the " effective cosmological parameters" . After that it turns out that so obtained equations 
of motion in general differ from the equations with the same parameters in the models based on the cosmological 
principle of uniform and isotropic Universe. If the difference looks like an effective dark energy contribution, then one 
can try to explain the accelerated expansion of Universe in frames of General Relativity without any dark energy. 

The present review considers all the three above listed possibilities. It aims to present a comprehensive description 
for both the kinematics of the cosmological expansion and the dynamical description of the process. We concentrate our 
attention on the acceleration of the cosmological expansion. It should be emphasized that the equivalence principle 
directly links the acceleration with the gravity nature and space-time geometry. Therefore its value is extremely 
important for verification of different cosmological models. It is the quantity which to great extent determines the 
ultimate fate of the Universe. 



This section is devoted to a particular way of Universe description, called the " cosmography" [l^ , entirely based 
on the cosmological principle. The latter states that the Universe is homogeneous and isotropic on the scales larger 
than hundreds megaparsecs, and it allows, among the diversity of models for description of the Universe, to select 
a specific class of uniform and isotropic ones. The most general space-time metrics consistent with the cosmological 
principle is the Friedniann-Robertson- Walker (FRW) one: 



Here a(t) is the scale factor and r is coordinate of a point, which does not take part in any motion except the global 
expansion of the Universe. In some sense the main problem of cosmology is to determine the dependence a(t). 

The cosmological principle helps to construct the metrics of the Universe and to make the first steps in interpretation 
of cosmological observations. Recall that kinematics means the part of mechanics that describes motion of bodies 
irrelative to forces causing it. In that sense the cosmography is nothing but kinematics of cosmological expansion. In 
order to construct the principal cosmological quantity — the time dependence of the scale factor a(t) — one have to 
have the equations of motion (the Einstein equations of general relativity) and an assumption on material composition 
of the Universe, allowing to construct the energy-momentum tensor. The efficiency of cosmography lies in the fact 



II. COSMOGRAPHY — THE KINEMATICS OF EXPANDING UNIVERSE 
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^ generally speaking, "epicycles" here stand not only for additional terms in the Hilbert-Einstcin action, but also for such new substances 
as dark energy and dark matter 
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that it allows to verify arbitrary cosmological models obeying the cosmological principle. Modifications of general 
relativity or introduction of new components (dark matter or dark energy) will certainly change the dependence a(t), 
but will not affect the relations between the kinematic characteristics. Velocity of the Universe expansion, determined 
by the Hubble parameter H{t) = a{t)/a{t), depends on time. The time dependence of the latter is measured by the 
deceleration parameter q{t). Let us define it by Taylor expansion of the time dependence of the scale factor a(t) in 
vicinity of the current moment of time to: 



a{t) = a {to) + a {to) [t - to] + ^'d {to) [t - to]^ + 
Let us rewrite the relation in the following form 



a{t)_ 
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where the deceleration parameter reads 
d{t)a{t) d{t) 1 
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As will be shown below, accelerated growth of the scale factor occurs at g < 0, while accelerated growth of the 
expansion velocity H > corresponds to q < —1. When the sign of the deceleration parameter was first defined, it 
seemed evident that gravitation — the only force governing the dynamics of Universe — dumps its expansion. A 
natural desire to deal with a positive parameter predetermined the choice of the sign. Afterwards it turned out that 
the choice made is inconsistent with the observed dynamics of expansion and it is rather an example of historical 
curiosity. 

For more detailed kinematical description of cosmological expansion it is useful to consider the extended parameter 
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Note that the latter four parameters are dimensionless. Derivatives of lower parameters can be expressed through 
higher ones. So, for example, 
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dln(l + z) 

Let us Taylor expand the time dependence of the scale factor using the above introduced parameters 
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In terms of the same parameters the Taylor series for redshift reads 
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Let us cite few useful relations for the deceleration parameter 
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For the case of single component liquid with density p one has 
,(a) = -l-^. 

The derivatives 4?, ^^iM--, ^zr and can be expressed through the parameters q and j solely 



(8) 



dH^ 

dz 
d'^H 
dz'^ 

11^ 

d^H 



1 + 9 



dz^ 



-H 



l + z 
J_ 

(1 + ^) 
H 

J)Tz) 

H 

OTzY- 



3 (3g2 + 3g3 _ 4gj - 3j - s) ; 



(-129^ - 24q^ - 15q^ + 32qj + 2bq^j + 7qs + 12 j - if + 8s + l) 



Let us present for convenience the useful relations for transformations from higher time derivatives to derivatives with 
respect to red shift: 
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Derivatives of the Hubble parameter squared with respect to the red shift i = 1, 2, 3, 4 expressed through the 

cosmographic parameters take the form 
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dz^ (1 + z)= 
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Time derivatives of the Hubble parameter can be expressed through above defined cosmographic parameters H, q, j, s, I 
as 

H = -H\l + q)- (f2) 

H = H^j+iq + 2)- 

H = H*[s-4.j -3q{q + A)-6]; 

H = H'^[l-5s + lQ(q + 2)j + 30{q + 2)q + 24]. 

From the relations one can easily sec that the accelerated growth of the acceleration velocity H > corresponds 
to the case g < — 1. 

The Hubble parameter, as can be seen from the relation ([5]), is linked to the deceleration parameter with the integral 
relation 



H = Hq exp 



[q{z') + l]dlii{l + z') 







It immediately follows that in order to calculate the principal characteristic of the expanding Universe H(z) one needs 
the information on dynamics of cosmological expansion, coded in the quantity q{z). 

III. DYNAMICS OF COSMOLOGICAL EXPANSION BRIEFLY 

Dynamics of Universe in frames of General Relativity is described by the Einstein equations 
1 

2' 



The energy- momentum tensor T^i, describes the space distribution of mass (energy), while the curvature tensor com- 
ponents i?^^ and its trace R are expressed through the metric tensor g^^ and its derivatives of first and second order. 
The Einstein equations generally are non-linear and hard to analyze. The problem is simplified when considering 
mass distribution with special symmetry properties provided by the metric. For homogeneous and isotropic Universe 
described by the FRW-metric, the Einstein equations reduce to the system of Fricdmann equations: 



2 

a ' 



SttG k 

= — P-^, (13) 

d AnG , „ ^ ,^ 

- = —{p + Sp). (14) 

a 6 

Here p and p are total energy density and pressure respectively for all the components present in the Universe at 
the considered time moment. 

The system (|13ll4p is sufficient for complete description of the Universe dynamics. The conservation equation 

p + 3H{p + p) = 0, (15) 

follows from the Lorentz-invariance of the energy-momentum tensor = and, as can be easily seen, it is nothing 
but the first principle of thermodynamics for ideal liquid with constant entropy dE -\- pdV = 0. We remark, that this 
equation can be derived from the Friedman equations (I13ll4p . 
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Let us now consider relative density for ith component of the Universe: 



2 



n.^f^, P..^, (i6) 

Setting the relative curvature density by definition flk = — qis^ ^ ; the first Friedmann equation can be presented in 
the form 

= 1. 

i 

In order to solve the Friedmann equations one has to define the matter composition of the Universe and construct 
the state equation for each component. In the simplest linear form the state equation reads 

Pi = Wip.^. (17) 

Solving the Friedmann equations for the case w = const, k = one obtains 

a{t) oc (t/to)^^, p oc a-^^^+'^l 

(The scale factor is normalized by the condition a (to) = !■) The above cited solutions exist only provided w 7^ — 1, 
the otherwise case will be considered separately. For the Universe dominated by radiation (relativistic gas of photons 
and neutrinos) one has w = 1/3, while for the matter dominated case w = 0. As a result of such state equations, for 
the matter dominated case one gets 



a{t) (X [t/tof'^^ , p (X a 

which is easily interpreted as the particle number conservation law. For the radiation-dominated case the solution 
reads 

a{t) (X (t/to)^^^ , pcca^'^. 

The latter result is a consequence of the fact that the energy density of radiation decreases as a^"^ due to growth 
of occupied volume (as Universe expands) and additionally as due to the redshift. Note that from the equation 
()15p it follows that p — const for w = —1. In the latter case the bubble parameter is constant leading to exponential 
growth of the scale factor 

a{t) cx e^*. 

Therefore in the case of traditional cosmological components in form of matter and radiation with w — and w = 1/3 
respectively, the Universe expansion can only decelerate, i.e. a < 0. 

Using the definition of deceleration parameter, one can see that for flat Universe filled by single component with 
state equation p = wp 

Q = i(l + 3u;). 

In generic case {k — (0, ±1), P — ^Pi, P — Pi'Wi) one gets 

i i 

9 = ^ + ^11^*^- (18) 



Using (|16p . the latter relation can be rewritten in the form 



Since the deceleration parameter g is a slowly varying quantity (e.g. q ~ 1/2 for matter-dominated case and q = —1 
in the Universe dominated by dark energy in form of cosmological constant), then the useful information is contained 
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in its time average value, which is very interesting to obtain without integration of the equations of motions for the 
scale factor. Let us see how it is possible [2T'|. Let us define the time average q on the interval [0, to] by the expression 

9(^o) = - / qit)dt. 

Substituting the definition of the deceleration parameter 
da 0? / 1 \ 



it is easy to see that 

to Ho 



q{to) = -1 + -4r (19) 



or, equivalently, 

io = (20) 

As was expected, the present age of universe scales as but the proportionality coefficient depends only on 

the average value of the deceleration parameter. It is worth to mention that this purely kinematic result does not 
dependent on curvature of space in the Universe, neither on number of component filling it, nor on the particular 
type of the gravity theory. Let us somewhat reformulate the obtained results for the average deceleration parameter. 
For single-component fiat Universe the Friedman equations can be presented in the following form 



SttGp = 3^; (21) 

• 2 

SttGp = -2--%r. (22) 



For the single component Universe with the state equation p = wp, w = const the equation for the scale factor 



1 + 3u; 



= 0, 



has an exact solution 
a{t) = ao 



+ Hot 



(23) 



Now it immediately follows that 

q^^—^ const, to^-T-^y 24 
2 3(1 + w) 

In particular, one has for the cases of radiation {w = 1/3) , q = I and non-relativistic matted (w = 0) , q = 1/2. From 
the other hand, accounting that in the considered case q — q, the latter relation can be rewritten as to — Hq^/ (1 + q) , 
which coincides with the relation (|20p . The resulting expression ([Ml) can be presented in the following form 

(25) 

1 + 9 

where T, iJ, q are age of Universe, Hubble parameter and average deceleration parameter respectively. Since q is of 
order of unity, then it immediately follows from ([25)1 that on any stage of Universe evolution the Hubble time ^^o"^ 
represents the characteristic time scale. 

Further model-independent dynamical restrictions on the Universe kinematics can be derived from the so-called 
energy conditions [22, 24-27]. Those conditions based on rather general physical principles, impose restrictions on 
the energy- momentum tensor components T^i/. Specifying a particular model for the medium^ those conditions can 
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be transformed into inequalities, limiting possible values of pressure and density in the medium. In the Friedmann 
model the medium is an ideal fluid, where 

T^^u ^ ip + p) Uf^Ui, ~ pg^^, (26) 

where is the fluid four-velocity, with total density p and pressure p given, respectively, by 

3 / a? k\ \ ( 'a k 



8.GV^ + ^j' ^ = -8^(2- + - + -). (27) 

The most common energy conditions (see, e.g., |22l l24| - |26j ) reduce to 

NEC => p + p>0; 

WEC ^ p>0 and p + p>0; 

SEC ^ p + 3p > 0, mid p + p>0; 

DEC ^ p>0 and - p<p< p, 

where NEC, WEC, SEC and DEC correspond, respectively, to the null, weak, strong and dominant energy conditions. 
Because these conditions do not require a specific equation of state of the matter in the Universe, they provide very 
simple and model-independent bounds on the behavior of the energy density and pressure. Therefore, the energy 
conditions are among of many approaches to understand the evolution of Universe. From Eqs. (j27p . one easily 
obtains that these energy conditions can be translated into the following set of dynamical constraints relating the 
scale factor 

NEC + ^ + A > 0; 

a a 

WEC 4 + 4 > 0; (28) 
a"^ 

SEC ^ - < 0; 



DEC 



> 0. 



For flat Universe the conditions (|28|) can be transformed to restrictions on the deceleration parameter q 
NEC ^ q>-l; 

SEC ^ q>0; (29) 
DEC ^ q<2, 

while WEC is always satisfied for arbitrary real a{t). 

The conditions (|29|) . considered separately, leave a principal possibility for both decelerated {q > 0) and accelerated 
(g < 0) expansion of Universe. The sense of the restrictions for NEC in is quite clear. As it follows from the 
second Friedmann equation, the condition of accelerated expansion of Universe reduce to the inequality p + 3p < 0, i.e. 
accelerated expansion of Universe is possible only in presence of components with high negative pressure p < —1/3/9. 
The energetic condition SEC excludes existence of such components. As a result q > in that case. At the same 
time the conditions NEC and DEC are compatible with the condition p < —l/3p and therefore regimes with q < 
are allowed here. In conclusion of the section let us pay attention to an interesting feature of the expanding Universe 
dynamics. According to Hubble law, the galaxies situated on the Hubble sphere, recede from us with light speed. 
Velocity of the Hubble sphere itself equals to time derivative of the Hubble radius Rh = c/H, 



. d f 1 



c a 



dt^^-^ = 'dt[H)--H^[-a-^)-^'^' + '^^- (^0) 

In the Universe with decelerated expansion (q > 0) the Hubble sphere has superluminal velocity c{q + 1) and it 
outruns those galaxies. Therefore the galaxies initially situated outside the Hubble sphere, gradually enter inside it. 
Any observer in arbitrary point of the Universe will see that the number of galaxies constantly increases. In the case 
of accelerated expanding Universe (q < 0) the Hubble sphere has subluminal c(l — |g|) and it lags behind the galaxies. 
Therefore now the galaxies initially contained in the Hubble sphere will gradually leave it and become unaccessible 
for observation. Should we consider those unobservable galaxies as physically real? The distinction between physics 
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and metaphysics is the possibihty of experimental verification for the physical theories. Physics have no deal with 
unobservable objects. However, the science constantly expands its boundaries, thus including to consideration more 
and more abstract concepts, once been metaphysical before: atoms, electromagnetic waves, black holes... The list can 
be continued further. 

We are most likely inhabitants of accelerated expanding Universe. Likewise in the decelerated expanding one, there 
are galaxies so distant from us that no signal from them can be detected by terrestrial observer. However, if the 
cosmic expansion accelerates, then we recede from those galaxies with superluminal velocity. Therefore, if their light 
have not reach us till now then it will never come. Such galaxies are unaccessible for observation not only temporarily, 
they are forever unobservable. Such "never observable galaxies" descent from the same "Big Bang" like our Milky 
Way. Should we attribute them to physical objects or to metaphysics? Those who consider the science fiction as a 
realization of most unbounded fantasy, are absolutely wrong. Compared to modern cosmology, the science fiction is 
dull and lacks imagination. 

Accelerated expansion of the Universe was first introduced in cosmological models with creation of the inflation 
theory. It was developed to fix multiple defeats of the Big Bang model. It turns out that in order to resolve most of 
them it is sufficient to have exponentially fast accelerated expansion of Universe at the very beginning of its evolution, 
during just about 10~^^s. The simplest way to obtain expansion of such type is to consider dynamics of Universe with 
scalar field. The inflation theory is formulated in many ways, starting from the models based on quantum gravity and 
high-temperature phase transitions theory with supercooling and exponential expansion in the false vacuum state. 
To illustrate the main ideas of the theory let us consider flat, homogeneous and isotropic Universe filled by the scalar 
field ip with the potential V{p)^ independent of the coordinates. The first Friedmann equation (|13p in the considered 
case takes the following form 

The conservation equation ()15p . written in terms of the scalar field, takes the form of Klein-Gordon equation, and in 
the case of non-stationary background one obtains: 

p + iHLp + V'{Lp) =Q. (31) 

In the fast expanding Universe the scalar field rolls down very slowly, like a ball in viscous medium, while the effective 
viscosity turns out to be proportional to the expansion velocity. In the slow-roll regime 

H(p:^(p, V{ip):s>(p'^. 

In the same limit the equations of motion take the form 

3HiP + V'{ip) =0; 



To be more specific, let us consider the simplest model of the scalar field with mass m and potential energy V{ip) 

1 



^Lp"^ . Right after the start of inflation ip <^ iHp] p^ <^ m?p'^^ therefore 



i—tp m^p = 0; 
a 

a 2jnip pK 



H . 

a Mpi V 3 

Due to fast growth of the scale factor and slow variation of the field because of strong friction, one has 



Mpi V 3 

For better clarity let us obtain the state equation for the scalar field in the slow-roll regime. For homogeneous scalar 
field in the potential V {ip) in the local Lorentz frame the non-zero components of the energy momentum tensor are 

Too = ^<y5^ + V{p)^ p^; T,j = (^p'^ ^ V {p?)j % = p^S^j. 



10 



In the slow- roll regime tp^ <C V{f) and, consequently, ~ —p^p. Therefore the energy-momentum tensor in the 
slow-roll regime approximately coincides with the vacuum one with p = —p. Using (jlSp and taking into account the 
fact that the inflation scenery with exponential growth of the scale factor leads to the flat space case, and therefore 
f2 = 1, one obtains that g = — 1 during the period of the inflation expansion. 

Therefore the scalar field provides rather natural way to obtain the accelerated expansion of Universe at least on 
early stages of its evolution. As the field intensity decrease in the slow-roll regime, viscosity falls down and the 
Universe comes out the inflation regime with exponential growth of the scale factor. 

We remark that the scalar field can provide decelerated expansion of Universe as well as the accelerated one. Near 
the minimum of the inflation potential the inflation conditions definitely violate and the Universe leaves the inflation 
regimes. The scalar field start to oscillate near the minimum. Assuming that period of the oscillations is much less 
than the cosmological time scale, one can neglect the expansion term in the equation (j3ip . and it easy to come up 
with effective state equation near the minimum of inflation potential. Let us rewrite the scalar field equation (l3ip in 
the form 

^ (^^) - + V'K = 0- 
After averaging over the oscillation period the first term turns to zero and then 

(^') - «) ■ 
The effective (averaged) state equation now reads 

{^V^) + (2U) ■ 

For the case of quadratic potential V (yi Lp^ one obtains ~ 0, which corresponds to the state equation of non- 
relativistic matter. 

We note also that the scalar field models are very wide-spread in modern cosmology, allowing to describe not only 
accelerated expansion, but much more complicated dynamics of Universe. Besides, most of the scalar field models are 
well motivated in the particle physics and concurring theories. 



IV. EVIDENCE FOR THE ACCELERATED EXPANSION OF UNIVERSE 



The Hubble law tells nothing about the magnitude, sign or very possibility of non-uniform expansion of Universe. 
The approximation where it works is insensible to acceleration. Investigation of the non-linear effects requires high 
red-shift data. If the observation detect deviations from the linear law, then the magnitude and sign of the deviation 
allow to judge about the sign of the cosmological acceleration. If the detected deviation lies towards larger distances 
at fixed redshift then the acceleration is positive. The distance is estimated by brightness of the source in assumption 
that the considered kind of sources represent the standard candles — an ensemble of object with practically equal 
luminosity. Therefore the observed brightness of the objects depends only on the distance to observer. The supernova 
outbursts of type la (exploding white dwarfs) present an example of such kind of objects. As the white dwarfs 
have very low mass dispersion, their luminosity is practically the same. An additional advantage is the huge power 
of (~ lO'^^iy) , released in explosion. Thanks to that, they can be detected on distances compatible to size of the 
observable Universe. 

Given the emitted light intensity L, or internal luminosity of the object, and having measured the light intensity 
F, which have reached us, or observed flow, one can calculate the distance to the object. The quantity defined this 
way is called the luminosity distance cLl 

In order to determine the acceleration of the Universe expansion one needs to express the luminosity distance in terms 
of the redshift for the registered radiation. Let E be internal (absolute) luminosity of some source. A terrestrial 
observer registers the photon flow F. Increasing of the photon wavelength, and therefore decrease of its energy, in the 
expanding Universe during its pass from the source to the observer results in effective (apparent) luminosity of the 
source L — E/a{t). The conservation law for the energy emitted on time interval dt and absorbed on the interval dto, 
implies 

FAnr'^dto = Edt = La{t)dt, (33) 
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where r is the comoving distance between the source and the observer at time moment to, which coincides with the 
physical distance in normahzation of the form a (ip) = 1- Since the comoving distance between source and observer 
does not change, the conformal time interval dr] ~ dt/a between two hght pulses is the same at the point of emission 
and at the point observation 

dt dto 
a{t) ao ' 

Therefore from ((33|) follows that 

La^t) 



F 



Comparing this expression with the definition of the luminosity distance (|32|) . one finds 



a{t) 



{l + z)r. 



(34) 



(35) 



The obtained result is physically clear. In the expanding Universe the registered flow decreases by factor (1 + z)^ : 
first, due to increase of the photon wave length (1 + z) times and, second, due to increase of the time interval for 
arrival of the fixed energy portion also {1 + z) times. 

Let us now determine the comoving distance to the currently observed light source as the function of its redshift. 
The equation of motion for photon is ds^ = 0. Let us consider a radial trajectory with the observer at the origin of 
the coordinates. In the case of spatially flat metrics one has 



ds^:^a^{t) {drj^- 
Accounting that 

^ dr] dt da ^ 
dt da dz 



dr^) = 0. 



dz 



H{zy 



one finds 

t{z) = 



dz' 



H{z') 

Therefore for spatially flat Universe 
dz' 



di = (1 + z) 
In general case. 



H{z'y 



(36) 



(37) 



(38) 



(39) 



dL{.z) = c(l + z)(l - UoY^'^H^^S 



.1/2 



dz' 



H{z') 



Ho ^ , (40) 



where 



S{x) 



sinx, r^o > 1; 
X, flo = 1; 
sinhx, fio < 1- 



The quantities Ha,fla = f^, (^f 



3Ho 
SttG 



1 + Z 

Ho 



dz' 



(41) 

refer to the present moment of time. For the multi-component flat case 

(42) 



lj:no.{i + z'f'+-^^ 

The relation (|40l) . can be rewritten in terms of the deceleration parameter. For the flat case it is 



dLiz) = il + z) 



(43) 
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Let us use normal definition for tfie absolute value fj, of the distance to a standard candle (say, SNela): 

^(z) = [mBiz) - Mb] = Slog (dL/Mpc) + 25. (44) 
Using one can link that quantity to the acceleration history q{z) via 

l + z 



H{z) = 25 + Slog 



HoMpc 



duexpy-J [1 + q{u)]dlnu ] . (4S) 



Here Mb and ttib are absolute and apparent stellar magnitude of the source respectively. The latter expression 
represents a fundamental relation, linking the deceleration parameter history with the SNela measurements. 

We remark that the relation (j43|) is based solely on the FRW-metrics. It means that the acceleration-deceleration 
dilemma can be resolved without any assumption about applicability of general relativity. However, supernova 
observations cannot help to determine current deceleration parameter immediately. In order to interpret the data, we 
must know H{z) or q{z), which in turns requires to specify the dynamical equations and material composition of the 
Universe. 

It is useful to cite the expression for the luminosity distance with accuracy up to z^ terms: 



no 



i+(^)^ + o(.^) 



(46) 



where in the flat case go = + 3t«i). 

i 

It follows from that for small z the luminosity distance is proportional to the redshift, with the proportionality 
coefficient equal to inverse value of the Hubble constant. For more distant cosmological objects the higher order 
corrections for the luminosity distance depend on the current value of the deceleration parameter go, ore, in other 
words, on quantity and type of components, filling the Universe. 

The expression for luminosity distance in the next order correction over the redshift reads 

dL{z) = -^[l + i(l-9o)2-i(l-go-3go'+jo)2'+ (47) 

+ ^ (2 - 2go - ISgg - Ihql + Sjo + 10jo9o + so) ^3 + 0{z% 

Let us briefiy discuss the methodology for verification of the cosmological models using the type SNela supernova 
outbursts. Let us consider more closely [1, Q two supernovae 1992P at low redshift z — 0.026 with m — 16.08 and 
1997ap at high redshift z — 0.83 with m = 24.32. As indicated earlier, dL{z) ~ z/Hq for z ^ 1. Using pi)) we find 
M = -19.09. Then the luminosity distance of 1997ap is obtained by substituting m = 24.32 and M = -19.09 for (gl]) 

HodL =i 1.16 for z = 0.83. (48) 

From the other hand, for the case of Universe filled by non-relativistic matter, using one finds 

HodL ~ 0.95. 



The latter result apparently contradicts the observations (l48t . 

Therefore, registration of outbursts for any kind of light sources with the same internal luminosity, i.e. standard 
candles, allows to determine the velocity of Universe expansion in different moments of its evolution. Comparing the 
obtained results with the predictions of different theoretical models, one can select the most adequate ones. Despite 
of the conceptual simplicity, the problem have been continuously facing many obstacles. Let us list just a few of 
them. The supernova outbursts are rare and random. To collect sufficient statistics one has to cover rather large part 
of the sky. The outbursts takes place for a limited time, therefore it is important to detect a supernova as soon as 
possible in order to follow the dynamics of its brightness variation. And of course the main problem is still connected 
to disputable applicability of the type la supernova for a standard candle. In the early ninetieth the United States 
created two groups for detection and analysis of the type SNela supernova outbursts: SuperNova Cosmology Project 
and High-Z SuperNova Search. It were their results jj, 01 which in the years 1998 and 1999 gave base for conclusion 
on accelerated expansion of Universe, so dramatically affected not only modern cosmology, but all physics as a whole. 
During the passed decade the results [1, 01 were multiple times reproduced with ever improving statistics. The main 
conclusion was always reconfirmed: comparable recently (at z ~ O.S) our Universe passed through the transition from 
decelerated expansion to accelerated one. Let us present in more details the analysis carried out by Riess [4^ . where 
he used the so-called golden set of SNela. The set contained 157 well studied SNela supernova with redshifts ranged 
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in 0.1 < z < 1.76. The analysis was based on the relation (|43l) which is exact expression for the luminosity distance 
in a geometrically flat Universe. In the case of linear two-parameter decomposition 

q{z) = qo + qiZ, (49) 

the integral in (^5)) can be exactly evaluated and the expression for the luminosity distance takes the form [4^ 

dL{z) = ^e'^'qf-'' [7 (91 " 9o, (1 + z)qi) - 7 (9i ^ 9o, gi)] , (50) 

where qo, qi are values of q{z), ^^^j^ at z = 0, 7 is the reduced gamma- function. Using ([50]) one can obtain information 
about qQ,qi and, therefore, about global behavior of q{z). The dynamical "phase transition" occurs at q{zt) = or 
equivalently at zt = —qo/qi- 

The second widely used parameterization is 

q{z)^qQ + qi— — . (51) 
1 + z 

Its advantage is in the fact that it is well-behaved at large z, while the linear approximation suffers of divergences. In 
such parameterization one has 

^^(^) = nJ^u e'^li^'"^'"^ [7 (91 + 90, 91) - 7 (91 + 90, 9i/(l + z))] . (52) 
(1 + z)Ho 

Now the parameter qi determines the correction to qo in the distant past: q{z) — qo + qi at z 3> 0. The likelihood for 
the parameters qo and qi can be determined from statistic, 

2,Tj \ sr^ {fJ-p.i{zi;Ho,qo,qi) - fj.Qi)'^ 

X (Ho,qo,qi) ^ 2^ 2 — \ — 2 ' (53) 

i l-i-a.i ^ V 

where ^ is the uncertainty in the individual distance modules, CTu is the dispersion in SNe redshifts due to peculiar 
velocities. The obtained results give evidence in favor of the Universe with recent acceleration (go < 0) and previous 
deceleration (qi > 0) with 99.2% and 99.8% likelihood. In the case of linear decomposition of the parameter q 
the transition from decelerated expansion in the past to current accelerated expansion took place at the redshift 
Zt = 0.46 ± 0.13. Unfortunately one should not take this result too serious: the linear approximation always leads to 
the transition, provided the two parameters have opposite signs. 

The consequent and statistically more reliable analysis qualitatively supported the above cited results. Of course, 
the quantitative results depend on the used set of supernova, but the main result remains the same: we live in 
accelerated expanding Universe, which passed from decelerated expansion to accelerated one in the recent past. 

Nowadays the method based on the supernova outbursts observations is undoubtedly leading. But it turned out to 
have competitors. A promising and completely independent alternative (and no way a surrogate) is the observation of 
angular diameter distances Da{z) for a given set of distant objects. The combination of the Sunyaev-Zeldovich effect 
with measurements of surface brightness in X-ray range provides the method for determination of angular diameter 
distances for galactic clusters H^]. The Sunyaev-Zeldovich stands for small perturbation of CMB spectrum produced 
by inverse Compton scattering of relict photons passed through clouds of hot electrons. Observing the temperature 
decrement of the CMB spectrum towards galaxy clusters together with the X-rays observations, it is possible to break 
the degeneracy between concentration and temperature thereby obtaining Da{z). Therefore, such distances are fully 
independent of the one given by the luminosity distance. 

For the case of spatially flat Universe described by FRW-metric, the angular diameter distance is 



, , , 1 r- du 1 

dA{z) 



{l + z)HoJo H{u) {l + z)Ho 



/ exp 


[-1 


10 





[{l + qiu')]d\n{l + u' 



du, (54) 



considered the 38 measurements of angular diameter distances from galaxy clusters as obtained through SZE/X- 
ray method by Bonamente and coworkers [47| Lima et al used a maximum likelihood determined by a statistics 
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where Daq ^ is the observational angular diameter distance, aoAoA^ the uncertainty in the individual distance, astat 
is the contribution of the statistical errors added in quadrature and the complete set of parameters is given by 
p = (iJo, <Zoj <Zi) • For the sake consistency, the Hubble parameter Hq has been fixed by its best it value Hq = 
80km/ sec- Mpc^^. In the case of linear parameterization the results are the following: best its to the free parameters 
are go — —1-35, qi = 4.2, zt = 0.32. Such results favor a Universe with recent acceleration (go < 0) and a previous 
decelerating stage {dq/dz > 0) . In the case of the parameterization 

q{z) go + gi 



1 + z' 



one gets go = —1.43, gi = 6.18, zt = 0.3. In both cases the results well agree with the ones obtained using the SNela 
data. 

It was recently shown that the so-called luminous red galaxies (LRG's) provide another possibility for direct mea- 
surements of the expansion velocity [4^ |4§, [sH] • Idea of the method is to reconstruct the Hubble parameter from the 
time derivative of the redshift 

H{z) = --^^. (56) 
1 + z dt 

The derivative can be found in measurements of the "age difference " between two passively evolving galaxies at 
different but close redshifts. The method was realized for 0.1 < z < 1.75. We stress that the range includes the 
interesting for us transition region, where z ^ 0.5. The results of analysis for available data (56| agree with the ones 
based on the supernova and angular diameter. In nearest future a set of some 2000 passively evolving galaxies is going 
to be measured in the range < z < 1.5. Those observations will give about 1000 values of H{z) with 15% accuracy, 
provided the age of galaxies is determined with 10% precision. 

Concluding the present section we present a novel and promising direction to research the history of cosmological 
expansion of Universe. We remind that astrophysicists initially used cepheids — the variable stars with intensity 
proportional to period of brightness variation — as the standard candles. A typical example of cepheids is the Polar 
Star, the brightest and closest to Earth variable star with period of 3.97 days. The cepheids are perfect standard 
candles for galactic distances. They allowed to determine size of our Galaxy and distance to our closest neighbor — 
the Andromeda galaxy. Research of the Universe dynamics involves principally larger scales, and thus it requires to 
use essentially more powerful standard candles. Recall, that the cosmological principle, which postulates homogeneity 
and isotropy of the Universe, and which all our above used equations of Universe dynamics are based on, is valid on 
scales of 100 Mpc and larger. Using considerably more powerful radiation sources as the standard candles, namely type 
la supernovae, already allowed us to advance remarkably deeper to the history of Universe. However the possibilities 
of the new standard candles are also limited. Up to present time the type la supernovae were observed only up to 
the redshift z < 2, though more reliable reconstruction of the cosmological expansion history requires even higher 
redshifts and therefore more powerful standard candles. It turns out that objects with required properties are already 
available! It is the so-called gamma-ray bursts(GRB) — giant energy emissions of explosive type with duration from 
tree to hundred second, observed in the hardest part of the electromagnetic spectrum. The typical energy, emitted 
by the gamma-ray burst is ~ 10^* erg which is one order of magnitude higher than at a supernova outburst. It is 
already compatible to the Sun rest mass! The events producing the gamma-ray bursts are so powerful that they can 
be observed even by naked eye, though they occur on distances of billions light years from the Earth! The energy 
yields in form of a coUimated flow, called a jet. Occurrence of the jets means that we can see only small number of 
all bursts that happen in Universe. Distribution of GRB durations has pronounced bimodal character. 

Origin of short GRB is probably linked to coalescence of two neutron stars, or a neutron star and a black hole. 
The events with longer duration are presumably due to creation of a black hole in collapse of massive stellar nucleus 
(/,25 Solar mass) with high angular momentum — it is the so-called collapsar model. The possibility to use the GRB 
as the standard candles is based on occurrence of the so-called "Amati relation", which relates the peak frequency of 
the burst with its total energy. This relation is a direct analog of the period-luminosity relation for cepheids. Strong 
dispersion (see FiglT|) yet limits applications of the GRB as the standard candles, but the possibility to advance to 
the region of considerably higher redshifts makes this direction potentially very attractive. 



V. SCALE FACTOR DYNAMICS IN SCM 



Now let us consider evolution of the deceleration parameter in Standard Cosmological model (SCM). Recall that the 
Big Bang model included only two substances - matter and radiation - which could give only decelerated expansion 
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FIG. 1: 15 GRB with measured redshifts and collimation angle(compared to la supernova data in SCM) 



of the Universe. According to SCM, the Universe is presently dominated by the dark energy — a component with 
negative pressure. It is the component responsible for the observed accelerated expansion of Universe. Let us determine 
the redshift and time moment of the transition to the accelerated expansion, i.e. find the inflection point on the curve 
(see Fig. [2]), describing the time dependence of the scale factor The second Friedmann equation in the SCM (ITil) can 
be transformed to the following 



d 1 2 
a = 2^° 



{1 + zY 



(57) 



where D,a_o and f2mo are present values of the relative density for dark energy in form of cosmological constant and for 
matter respectively. Now it follows for the redshift at the point of transition from decelerated expansion to accelerated 



1/3 



- 1. 



(58) 



With the SCM parameters JIao ~ 0.73, ilmo — 0.27 one gets z* ~ 0.745. Note that the result ([58|) can be obtained 
using the fact that for a multi-component Universe with state equations pi = WiPi respectively, the deceleration 
parameter (fT8|) equals (recall that 51 = ^ il^ = 1 in the flat case ) 



1 



AO 



(59) 



The condition q = allows to reproduce the relation (l58t . Let us discuss its asymptotes. For early Universe {z — oo), 
filled by components with positive pressure , q{z ~> oo) = i, i.e. as it was expected the expansion is decelerated, 
while in the distant future with domination of the cosmological constant the accelerated expansion is observed: 
g(z— 1) = —1. The latter result is a trivial consequence of the exponential expansion a oc e^* in the case of 
Universe dominated by dark energy in the form of cosmological constant. 

We remark that the dependence q{t) can be immediately obtained from the definition q = —d/aH^, using the SCM 
solutions for the scale factor: 



a{t) = 
A = 



f^AO 



^A 



g-^O "AO • 



(60) 
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FIG. 2: Time dependence of the scale factor in SCM. 



It results in the following 



1 - 3 tanh^ ( — 



(61) 



The dependence q{t) is presented on FiglH Note that the time asymptotes at t — > and t — > oo correspond to the 
above obtained ones at z — cx) and z — > — 1. 

Let us now determine the time of the transition to the accelerated expansion. Inserting the relation (1601) . one 
obtains 



-2,-1/ 



t{a) = ^n^^'^H„ ^arsh 



\ ^mO ) 



(62) 



Transforming from the redshift to the scale factor a* = (1 + 2*) ^ = 
t* = (a*) = ^n^y^H-^ sinh^i (1/2) ~ 5.25 Gyr. 



1/3 



one gets 



(63) 



Because of physical importance of the obtained result we present another one, probably the simplest, interpretation 
of it. If the quantity aH = a growth, then a > 0, which corresponds to accelerated expansion of Universe. According 
to the first Friedmann equation. 



aH 



0.73a3 + 0.27 



It is easy to show that the function under radical starts to grow at a* ~ 0.573, which corresponds to z* — 0.745. It 
is interesting to note that the transition to accelerated expansion of Universe {z ~ 0.75) occurred remarkably earlier 
than the dark energy started to dominate {z ~ 0.4). 
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FIG. 3: Time dependence of the deceleration parameter q{t) in SCM. 



In the SCM the dark energy exist in form of cosmological constant with the state equation pA = —pA, i-e. the 
state parameter equals w\ = —1. A natural question rises; what is the limiting value of that parameter, which still 
provides accelerated expansion of Universe in the present time? As we have seen above, the condition for accelerated 
expansion is the following ^ (pj + 3pi) < 0. In SCM it is transformed to the following 

i 

Woe < -^^Bb; ^"de < -0-46. 

Of course, a substance with such state equation is not the cosmological constant, and it can be realized, for instance, 
with the help of scaler fields (see next section). 

Now let us estimate absolute magnitude of the cosmological acceleration using the SCM parameters. Taking time 
derivative of the Hubble law, one gets 



V 



R. 



(64) 



The time derivative of the Hubble parameter is 



a 



Therefore 



H + H^ = - = —{Pm- 2pa) = — [p^- 2'^'" ) = ^ ( "a - 

And finally we get the analog of the Hubble law for the acceleration V 



V = HR; H = H^ {nA--fl 



(65) 



In the present time (0^ = flmO, = ^ao), and on the distance, for example, R = 1 Mpc one gets 
V ~ 10-"cm sec"^. 



(66) 
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A possible way to observe such effect is based on the fact that the redshift for any cosmological object slowly changes 
due to the acceleration (or deceleration) of the Universe expansion. Let us estimate the magnitude of that effect. By 
definition 

a(io) _ 
a{t) 

where t is the time of emission and to is the detection time, and it follows that 

(67) 



dz _ a{to)a{t) - a{to)d{t)^^ 
dto a?{t) 



Taking into account that ^ = ^^K^ — jp^i^i one gets 
. _ d{to) a{to)d{t) 1 

a{t) a^t) l + z' ^ ' 

Thus, the rate of redshift variation for the light emitted at time t and registered at the present time is determined 
by the relation 

dz 

z=^=Ho{l + z)-H{t). 
dto 

In the SCM 

H = Ho{nmo{l + zf + Oao)'/^. 
Therefore for variation of the redshift on the time interval At one obtains 



Az = zAt = Ho 



l + z-(n„o(l + 2) +^Ao)'^' • (69) 



Note that for the two limiting cases rt\o = 1, ^mo = (accelerated expansion) and ^ao = 0, flmo = 1 (decelerated 
expansion) , as was expected Az has opposite sign. For the SCM parameters the redshift variation Az and the velocity 
increment AV for a source with the redshift z = 4: and duration of observation Ato = 10 equal respectively 

Az 

A^pslO"^ AV = c « 6 cm/ sec. 

l + z ' 

The result discourages by its smallness. However, taking into account fast progress in precision of observable cosmology, 
we should not despair. Let us give an example. Today the list of exoplanets - planets out the Solar system - counts 
more than 300 items. The most successful method to detect the exoplanets is to measure radial velocity of the stars. A 
star accompanied by an exoplanet experiences velocity oscillations " to us and back from us" , which can be measured 
observing the Doppler shift of the star spectrum. For the first look it is impossible: affected by the Earth, velocity 
of the sun annually oscillates with amplitude of centimeters per second. Even effect of Jupiter produces variations 
of only meters per second, while thermal broadening of the spectral lines for the star corresponds to dispersion of 
velocity of order of thousands kilometers per second. It means that even in the case of Jupiter one have to measure 
shift of the spectral lines on a thousandth of its width. It seems incredible, but the task was successfully worked out. 



VI. DYNAMICAL FORMS OF DARK ENERGY AND EVOLUTION OF UNIVERSE 

The cosmological constant represents one of many possible realizations of the hypothetic substance called the dark 
energy, introduced for explanation of the accelerated expansion of Universe. As we have seen above, such substance 
have to have the parameter w in the state equation p = wp, which satisfies the condition w < —1/3 (in absence of 
other components). Unfortunately the nature of dark energy is absolutely unknown, which produces huge number 
of hypotheses and candidates for the role of fundamental contributor to the energy budget of the Universe. We told 
many times about impressive progress of observational cosmology in the last decade. However we still cannot answer 
the question about time evolution of the state parameter w. If it changes with time, we must seek an alternative to the 
cosmological constant. For quite a short time a plenty of such alternative (with respect to w = —1) possibilities were 
investigated. The scalar fields, formed the post-inflation Universe, are considered to be one of the principal candidates 
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for the role of dark energy. The most popularity acquired the version of the scalar field with appropriately chosen 
potential V (if) . In such models, unlike the cosmological constant one, the scalar field is truly dynamical variable, 
and dark energy density depends on time. The models differ by the choice of the scalar field Lagrangian. Let us start 
from probably the simplest dark energy model of such type, named the quintessence. We define quintessence as the 
scalar field ip in the potential V {ip) , minimally coupled to gravity, i.e. experiencing effect of the space-time curvature. 
Besides that we chose canonical form for the kinetic energy term. The action for such field takes the form 



2 dxf^ dxjj 



(70) 



where g = dct g^^ . Equation of motion for the scalar field is obtained by variation of the action with respect to the 
field, 



1 



--da 



'99 



dip 
dxu 



dp 



(71) 



In the case of flat Friedmannian Universe, namely in FRW- metric ([T]) one obtain for homogeneous field pit) the 
following 



-m^^v'{ip) = ^, 



(72) 



where y'{p>) = This relation is sometimes called the Klein-Gordon equation. 

The energy- momentum tensor for the scalar field can be found by variation of (1701) with respect to metrics g^^ 
resulting in the following 



8S dp d(p 



g Sgf^" dx^ dxu 



9tJ.uL. 



(73) 



In the case of homogeneous iieldip{t) in locally Lorentzian reference frame, where metricsg^j/ can be replaced by the 
Minkowski one, we obtain the density and pressure of the scalar field in the following form 



= Too = ^^'^ + V{(p): 



= = ^p'^ -V{p). 



(74) 



The Friedmann equations for the flat Universe, filled by the scalar field, take the form 



inG 



ip' + V{cp) 



(75a) 



H = -AnGip'. (75b) 

The system (|75al - [75b|) should be complemented by the Klein-Gordon equation (|72|) for the scalar field. Note, that 
the latter can be obtained from the conservation equation for the scalar field 

p^ + 3H{p^+p^)^0. (76) 

by substitution of the expressions (1741) for energy density and pressure. 
Using (TM]) . one get the state equation for the scalar field 



_p^- 2V 



1 ,^2 



(77) 



As one can see, the parameter of the state equation for the scalar field Wip varies in the range 

-l<w^ < 1. (78) 
The state equation for the scalar field is conveniently transformed into the following 
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The function w{x) increases monotonously from the minimum value Wmin = — 1 at x = up to the maximum 
asymptotic value Wmax = 1 at a; — ^ oo, corresponding to = 0. In the slow-roll limit, x ^ 1 (((P' ^ V{lp)) the scalar 
field behaves as the cosmological constant, — —1. It is easy to see that in such case p^p = const. The opposite 
limit X 1 (</j ^ y ((/?)) correspond to the hard matter — 1. In that case the scalar field energy density evolves 
as Pip (X a~^. The intermediate case x ~ 1, p ~ corresponds to non-relativistic matter. 
Transforming ()76p into the integral form 



poexp 



-3 / (l + wp)^ 



(80) 



one finds that in the general case the scalar field energy density behaves as 

P0cxa~", 0<TO<6. (81) 

The value Wtp = — 1/3 represents boundary between the regimes of decelerated and accelerated Universe expansions. 
Thus, the case of accelerated expansion takes place for < m < 2. A natural question arises: what are the potentials 
for scalar fields that can provide the accelerated expansion of Universe? The question can be posed another way: 
what potentials can produce the quintessence applicable for the role of dark energy? Consider simplified version of 
the problem. Let us find the scalar field potential resulting in power law for the scale factor growth: 

a{t) oc tP. (82) 

For accelerated expansion the condition p > 1 must be satisfied. Recall that p = 2/3 for the case of non-relativistic 
matter dominated Universe, and p = 1/2 for radiation-dominated one, so in both cases the expansion will be decel- 
erated. Making use of the Friedmann equations, one can express the potential V{(p) and (p in terms of H and H. It 
allows to write down the system of equations, which describe the parametric dependence V on if : 





(83b) 

Excluding the time variable with the help of the relation — -y/^lni, for the case of the power law one finds 



Vi^)^Voe.J-J'-^^). (84) 
V V P rnpij 

The obtained result implies that the scalar field in the potential (IMl) under condition p > 1 can be treated as the dark 
energy, i.e. it provides accelerated expansion of Universe. As we can see above, in the quintessence model the required 
dynamical behavior is achieved by the choice of the scalar field potential. Another dark energy model, realized by the 
scalar field with modified kinetic term, is called the fc-essence [13, HH . 

Let us define the quantity X = ^9^" and consider the action for the scalar field of the form 



S = 



j d-^Xy/^L{ip,X). (85) 



where L is generally speaking an arbitrary function of the variables tp and X. The traditional scalar field action 
corresponds to the case 

L{ip,X)=X-V{p). (86) 

We restrict our consideration to the following subset of Lagrangians: 

L{p,X)=K{X)-V{^), (87) 

where K{X) is positive defined function of the kinetic energy X. For description of homogeneous and uniform Universe 
we should choose X — i<y9^. Using the standard definition ([7i]), one finds 

Pp^Liip,X) = K{X)-V{cp); (88a) 
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p^^2X^^^~KiX) + Viip). (88b) 
Accordingly, the state equation for fc-essence takes the form 

2X^-K{X) + VM 

We demonstrate the main features of fc-essence in example of the simplified model [s^l , where the Lagrangian has the 
form L = F{X). Such model is called purely kinetic fc-essence. In that case one has 

dF 

p^ = 2XFx'F- Fx = -gj^; (90a) 
F; (90b) 



The equations of motion for the field can be obtained either from the Euler-Lagrange equation for the action ((85)) . or 
by substitution of density and pressure expressions ()90al and I90b|l respectively into the conservation equation for the 
fc-essence. It results in the following 

Fx'P + Fxx^^^ + 3HFx^ = 0, (91) 
or in terms of the kinetic energy X 

(Fx + 2FxxX) X + 6HFxX = 0. (92) 
The latter equation can be exactly solved 

XF^ = fca~^. (93) 
where constant fc > 0. 

The solution (|93p X{a) has one important property: all principal characteristics of fc-essence {p^,p,p, w^; see (|90al - 
I90cp ) as functions of the scale factor are completely determined by the function F{X) and do not depend on evolution 
of other energy densities. All the dependence of fc-essence on other components appears only by means of a{t). But this 
dependence is trivially a(t) cx ptot, and it is determined by the dominant component in the energy density. Solutions 
of such type are called the tracker solutions and its occurrence allows to approach the solution of the coincidence 
problem. It can be shown [ssj . that such property is present not only in purely kinetic fc-essence, but also its general 
case. 

Extensive set of available cosmological observations shows that the parameter w in the state equation for dark 
matter lies in narrow vicinity of w = —1. We considered above the region — 1 < w < —1/3. The lower bound of the 
region w — —1 corresponds to the cosmological constant, and all remaining interval can be realized by scalar fields 
with canonic Lagrangian. Recall that the upper bound w — —1/3 appears because of requirement to provide the 
observed accelerated expansion of Universe. But can we go beyond the bounds of that interval? 

It is rather difficult question for the component of energy we know so little about. General relativity usually 
impose certain restrictions on possible values of the energy- momentum tensor components, which are called the 
energy conditions (see section HIH) . One of the simplest restrictions of such type is the condition p + p > 0. Its physical 
motivation is to prevent the instability of vacuum. Applied to dynamics of Universe this condition requires that 
density of any allowed energy component should not grow as Universe expands. The cosmological constant with 
PA = represents the limiting. Taking into account our ignorance about the nature of dark energy, it is reasonable 
to pose a question: may such a mysterious substance differ so much from already known "good" energy sources and 
violate the condition p -I- p > 0? Accounting the requirement of positive defined energy density, valid even for dark 
one (it is necessary to make the Universe flat), and negative pressure (in order to explain the accelerated Universe 
expansion), the above mentioned violation will result in w < —1. Some time ago such component, called the fantom 



22 



energy, attracted attention of physicists '86']. Action for the fantom field minimally coupled to gravity, differs only 
by sign of the kinetic term from the canonic action for the scalar field. In that case the energy density and pressure 
of the fantom field are determined by the expressions 

P^ = roo = -i^2 + y(^); p^=Tu = ~\ip^-V{ip), (94) 
and state equation takes the form 

If (f^ < 2V{lp), then < —1. As an example let us consider the case of Universe, containing only non-relativistic 
matter (w — 0) and phantom field (w^ < — 1). Densities of the two components evolve independently: pm oc a^^ and 
Pip oc a^^'^^^™*'-'. If matter domination ends at the time moment then the solution for the scale factor at t > tm 
reads 



a{t) = a{Un) 



+ (l+Wp) 



(96) 



It immediately follows that for w,p < — 1 at the time moment tgn = (i+w ) ^™- scale factor and whole set of the 
cosmological characteristics (such as scalar curvature and energy density of the fantom field) for the Universe diverge. 
This catastrophe was called the Big Rip, which is precede by a specific regime which is called super-acceleration. Let 
us explain the origin of the super-acceleration regime on a simple example. Consider differential equation 

In the case A > the equation (1^51) realizes the non-linear positive feedback. Fast growth of the function a;(t) leads 
to the "Big Rip" (divergence of the function to infinity) on finite time period. Indeed, the general solution of the 
equation reads 

where B is the integration constant. At t = —B the Big Rip occurs. 

It is easy to see that the model (|97l) represents a particular version of the Friedmann equation for w^p < — 1. Since 
p,p cx a~3(i-i-«)^)^ then the first Friedmann equation can be presented in the form 



a 



= Aa-i(i+"'^)+i. (99) 



For instance, with Wip — the equation (l99l) exactly coincides with ((98)) 



VII. /(i?)-GRAVITY, BRANEWORLD COSMOLOGY AND MOND 

Although SCM can explain the present accelerated expansion of Universe and well agrees with current observational 
data, theoretical motivation for this model can be regarded as quite poor. Consequently, there have been several 
attempts to propose dynamical alternatives to dark energy. We considered them in the previous section. Unfortunately, 
none of these attempts are problem- free. An alternative and more radical approach is based on assumption that there is 
no dark energy, and the acceleration is generated because of gravity weakened on very large scales due to modification 
of general relativity. In frames of such broad approach three main directions can be selected: f{R) gravity, braneworld 
cosmology and modified Newtonian dynamics (MOND). Let us briefly consider those alternatives to SCM from the 
point of view under our interest: slowdown or speedup? 



A. f{R)- gravity 



Theory of f{R) gravity is created by direct generalization of the Einstein-Hilbert action with the replacement 
R — > f{R)- The transformed action reads 



23 



For the generalization we take the function f{R) depending only on Ricci scalar R, not including other invariants, 
such as Rfj^yR^^'^ . Motivation for such choice is the following: the action with f{R) is sufficiently generic in order to 
reflect main features of the gravity, and with all that it is simple enough to avoid technical difficulties in calculations. 
We remark that the function f{R) must satisfy the stability conditions 



/'(i?)>0, /"(i?)>0. 



(101) 



were the primes denote differentiation with respect to the scalar Ricci curvature R. Complete action for f{R) gravity 
thus reads 



2k 



d^Xy^fiR) + Srn (ffpi., "0) 



(102) 



Here "0 is common symbol for all material fields. Variation with respect to metrics after some manipulations gives 



f'{R)R^, - \f{R)g^u 



(103) 



where 



'2 dSn 



(104) 



and is the covariant derivative. 



associated with the Levi-Civita connection of the metric and □ = V^V. 



Leaving aside the complications connecting with variation procedure, let us concentrate our attention on the field 
equations (|103p . They represent partial differential equation of fourth order in metrics, as the Ricci scalar R already 
contains second order derivatives of the latter. For the action linear on i?, the fourth-order derivatives (the two latter 
terms on the left hand side of (|103|) ') turn to zero and theory reduces to the standard general relativity. 

Note that the trace of (I103P gives 



J'{R)R - 2f{R) + 3Df{R) = SttGT, 



(105) 



where T = g^^'^T^j^i, links R with T differential relation, unlike algebraic one in general relativity, where R = SttGT. 
This a direct hint on the fact that the field equation of the /(i?)-theory allow much wider set of solutions compared 
to general relativity. To illustrate such statement we remark that the Jebsen-Birkhoff's theorem, stating that the 
Schwarzschild solution represents the unique spherically symmetric vacuum solution, does not hold in the /(i?)-theory. 
Leaving aside the details, we note, that now T = does not require that i? = or even is constant. 

The equation (|105|) appears very useful for investigation of different aspects of /(i?)-gravity, especially for stability 
of solutions and weak field limit. In particular it is conveniently used for analysis of the so-called maximally symmetric 
solutions. It is the solutions with R = const. For R = const and T^i/ — the equation (|105p reduces to 



/'(i?)i? - 2/(i?) = 0. 



(106) 



For a given function f{R) this equation is algebraic for i?. If i? = is a root of the equation, then (|103p reduces 
in that case to Rf^^, = and maximally symmetric solution represents the Minkovski space-time. If otherwise the 
root of equation (jl06p equals R — const ~ C, then (jl03p reduces to — ^g^u and maximally symmetric solution 
corresponds to the space of de Sitter or anti-de Sitter (the cosmological constant in general relativity) depending on 
sign of C. 

Now let us consider how dynamics of Universe are immediately described in f{R) cosmology. Inserting FRW-metric 
into (I103P and using the stress-energy tensor in the form 



T^j^u ^ (p + p) Uf_,u„ ~ pgfj^^, 
one obtains 



(107) 



P+\{Rf' - f)-2.HRf" 



(108) 



2H + iH^ 



SttG 



p + R^f" + 2HRf" + Rf" + - (/ - Rf) 



(109) 
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As was mentioned above, the main motivation for /(i?)-gravity is in the fact that it leads to accelerated expansion of 
Universe without any dark energy. The easiest way to see it is to introduce the effective energy density and effective 
pressure 

Rf - f iHRf" 

Peff = J, — ; (110) 

Vef! = J, . (Ill) 

In the spatially flat Universe the density p^sf must be non-negative, as it follows (llOSp in the limit p — > 0. Then the 
equations (jlOSp , (|109|) take the form of standard Friedmann equations 

a' = ^P.//; (112) 



- = 3- iPeff + SPeff) ■ (113) 

The effective parameter Wef / in the state equation in that case equals 
Peff R'r + 2HRf" + Rf" + - Rf) 

Weff = — TTT, — ; ■■ • (tt4 

The denominator in (|114p is definitely positive, so sign of w^ff is determined by the numerator. In general case of 
metric /(_R)-model in order to reproduce (mimic) the equation of state for de Sitter (cosmological constant) with 
W(.f / = — I , the following condition must be satisfied 

/'" RH~R 

Let us consider two examples (regardless its realizability) . The firs is the function of the form f{R) cx i?". We can easily 
calculate Weff as function of n, if we assume the power law for time dependence of the scale factor a{t) = ao(t/to) 
(arbitrary dependence a{t) results in time dependence for Weff)- The result for n ^ 1 is 

Qin? — 9n + 3 ' 
and a in terms of n reads 
-2n2 + 3n - I 



(117) 



2 

Appropriate choice of n leads to desired value for Weff- For instance, n — 2 leads to 

We// = — 1, a = 00. (118) 

This result is expected from consideration of quadratic corrections to Einstein-Hilbert action, which were used in the 
inflation scenery by Starobinski j88| . 
The second example is 

2(n+l) 

In that case, assuming again the power law for time dependence of the scale factor, one gets 

"■"-' + 3(2nTlMnVl) - <™' 

The condition for accelerated expansion Weff < —1/3 for the two above considered models is transformed into the 
form 

=F n - 1 > 0, 

where the signs =p respond to the first and the second examples respectively. In particular, in the second case, for 
n = 1 we find that lUg// = —2/3. We remark that in such case the positive values of n imply the presence of terms 
inversely proportional to R, unlike the above considered case. 
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B. Braneworld Cosmology 



All models explaining the accelerated expansion of Universe share one common property: they involve one or 
another way to weaken the gravity. It is the negative pressure in SCM and in the scalar field models, transformation 
of the universal gravity law in the modified gravitation model or voids in the models with inhomogeneities. An original 
approach to suppress the gravity is realized in the Braneworld scenario [s^. According to latter we live on a three- 
dimensional brane, which is embedded in the volume of higher dimension (in the simplest case it is four-dimensional). 
All material fields are restricted to the brane, while the gravitation penetrates both the brane and whole embedding 
volume (see Lecture notes by R.Maartens [9J|). Higher dimension of the space where the gravity acts, results in its 
suppression {F (x R~'^^~^\ where D is the dimension of the space). 

According such scenario the equation of motion for the scalar field on the brane reads 

if + 3Hip + V'{ip) = 0. (121) 

where 

H^^^p(l + ^)+^+'; (122a) 

3Mpz^ V 2a/ 3 a4' ' 

p=l^^ + V{^). (122b) 
Here e is the integration constant, which transforms the volume gravity to the brane restricted one. The brane tension 

(5) 

a provides connection between the for-dimensional Planck mass Mpi and five-dimensional one Mp^ : 



r(5)3 
PI 

47r I a/ct 



Mp, = wAM^|. (123) 



The tension a also links the four-dimensional cosmological constant A4 on the brane with the five-dimensional volume 
one Af, by: 

47r / . 47r , , , , 

Remark, that the equations (I122al - [l22bp contain an additional term c>c p^/cr, whose presence is connected with the 
conditions on the brane- volume boundary. The latter term is responsible for dynamical enhancement of the decay 
experienced by the scalar field while it rolls down the potential. In such a case infiation regime can be realized even 
for those potentials which are too steep to satisfy the slow-roll infiation conditions in the standard approach. Indeed, 
the slow-roll parameters in the braneworld models in the limit V/a 3> 1 read: 

e~4£Fw^-^^ ; ri~2riFRw(^^ ■ (125) 

It follows now that the slow-roll infiation condition (e, ij <^ 1) is easier satisfied than in the FRW cosmology under the 
condition V/a S> 1. The infiation thus can take place for very steep quintessence potentials, such a.s V oc e~'^'^, V cx 
(p'" and others. This in turn let us hope that inflation and quintessence can be generated by the same scalar field. 
A radically different approach to provide the accelerated expansion of Universe was proposed in papers [9^ , [9l| , 

M- 

The sharp distinction of the DDG model from the RS one is in the fact that both the volume cosmological constant 
and brane tension are set to zero, while a curvature term is introduced into the brane action. The theory is based on 
the following action 

S = M'^f f R + Mpi^ f R+ I L^atter- (126) 
J bulk J hrane J hrane 

The physical sense of the term /f,^^„g R is unclear. Probably, it accounts for quantum effects caused by material 
fields, which lead to the same term in the Einstein action, as was first noted by A.D. Sakharov in its induced gravity 



theory [93 
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The Hubble parameter in the DDG braneworld takes the form 



^=J^ + ^ + f, (127) 



where 1^ = Mpi'^/u'-p^ is a new length scale, defined by the 4D Planck mass Mpi and 5D one Mpi. An important 
property of such a model is that the accelerated expansion of Universe does require presence of dark energy. Instead, 
since gravity becomes five dimensional on length scales R > Ic = 2H~^{1 — fi^)"^, one finds that the expansion of 
the Universe is modified during late times instead of early times as in RS model. 

More general class of braneworld models, which includes both RS and DDG models as particular cases, is described 
by the action 

S = M'^^f I {R- 2Ah) + I (Mpi^R - 2(7)+ [ L^atter- (128) 

J bulk Jhrane Jbrane 

For (T = Ab = the action (|128p transforms into the DDG model action (|126p . and for to = it reduces to RS-model. 

As was shown in [92], the action (|128p describes the Universe, where the accelerated expansion era occurs on later 
evolution stages with the Hubble parameter 



where 



n„i{i + zf + n„ + 2ni t n.,n{i + zf + vi„ + ni + sIa,, (129) 



Signs =F correspond to two possible ways to embed a brane in to the volume. As it was in DDG model, 1^ - if 

(5) 

Mp; ~ lOOMeV. On shorter length scales r <ti h and on early times we recover the general relativity, while on large 

length scales r ^ 1^ and long time periods the brane effects start to be important. Indeed, setting Mp^ — (fi; = 0) 
the equation (|129p reduces to the KCDM model: 

^^^^nr,,{l + zf (131) 

while for the case cr = A;, = the relation (|129p reproduces the DDG model. An important feature of the action (|128p 
is the fact that it generates the effective state equation Wef / < — 1 • It can be easily seen ^] from the expression for 
current value of the effective parameter in the state equation 



2go - 1 ^ Q,^ / ^ 

3(l-a^) 1 - r!„ Y f]„ + + + 17a, ■ ^ ' 

Taking the lower sign, one gets wq < —1. 

It is worth noting that in the latter model the accelerated Universe expansion phase is the transient phenomenon, 
which comes to end when the Universe returns to the matter dominated phase. 



C. MOND 



The so-called modified Newtonian dynamics (MOND) is sometimes considered as an alternative to the DM (dark 
matter) concept. MOND represents such modification of Newton physics, which allows to explain the flat rotational 
curves for galaxies without attributing to any assumptions on DM. MOND assumes that the second Newton law 
F — ma must be modified for sufficiently low accelerations (a <C ao) so that 

F = mafi{a/ao). (133) 

where fi{x) = x if x <C 1, and fi{x) = 1 if a; ^ 1. It is easy to see that such law leads to the modification of 
the traditional formula for the gravitational acceleration F — mgN {qn — GM/r^). Relation between the "correct" 
acceleration and traditional Newtonian one reads 



(134) 
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For a rotating point mass one gets a = /r (this purely kinematic relation is independent on choice of the dynamic 
model). But here a stands already for "correct" acceleration. It follows that 

v"" = GMao- (135) 

i.e. for sufficiently small accelerations the rotational curves for isolated body of mass M do not depend on radial 
distance r, where the velocity is measured. In other words, the considered theory predicts not only flat rotational 
curves but also the fact that individual halo, associated with a galaxy, has infinite extension. This prediction may cause 
a serious problem for MOND, as recent observations of galaxy- galaxy lensing support the result that the maximum 
halo extension is about 0.5 Mpc. The value uq, needed to explain the observations 

ao^lQ-^cm/s^. (136) 

is of the same order as cHo . It supports the hypothesis that MOND "can reflect the influence of cosmology on local 
particle dynamics" . 

Although the results of MOND agree very well with the observations of individual galaxies, it is unclear whether they 
will be as well successful for description of cluster structure, where strong gravitational lensing points out considerably 
denser mass concentration in the cluster center than was predicted by MOND. Another difficulty connected with 
MOND, is the fact that it is quite problematic to incorporate it in more general relativistic gravity theory. Presently 
it is not clear what predictions are given by the MOND-like theories for complicated gravity effects like the gravitational 
lensing. 



VIII. DYNAMICS OF UNIVERSE WITH INTERACTION IN THE DARK SECTOR 



SCM considers dark matter and dark energy as independent components of energy budget of Universe. The dark 
energy postulated in form of cosmological constant, introduced as early as by Albert Einstein in order to make possible 
the creation of stationary Universe model. The assumed absence of interaction between the two components means 
that the energy densities for each component obey independent conservation equations 

p^ + 3H {p, + p,) = 0. 

Coupling between the two components leads to modified evolution of Universe. In particular, the energy density 
for non- relativistic component will not evolve according the law a^"^, and the dark energy density (in form of the 
cosmological constant) will not remain constant any more. From one hand, such modification of the theory opens 
for us new possibilities for solution of principal problems in cosmology. So, for example, solution of the coincidence 
problem'^ reduces to appropriate choice of interaction parameters, which can provide satisfaction of the condition 

^ = 0(1), 

in the present time consistently with the condition a > 0. From the other hand, introduction of the interaction will 
modify the relations between the observable parameters. In particular, the modification will affect the fundamental 
relation between the photometric distance and the redshift, which the evidence of accelerated expansion of Universe 
is mainly based on. It imposes strict limitations both on the form of interaction and its parameters. 

Taking into account that dark matter and dark energy represent dominant components in the Unive rse, then from 
point of view of field theory it is naturally to consider interaction between them J67,, 169. - il74 , Il82l Il98j | . A ppropriat e 
interaction between the dark matter and dark energy can facilitate solution of the coincidence problem [I75l4l77 | 



Non-minimal coupling in the dark sector can considerably affect the history of cosmological expansion of Universe and 
evolution of density fluctuations, thus modifying the rate of the cosmological structure growt h. The possibi l ity of dark 
matter and dark energy to interact with each other is widely discussed in recent literature |l67l . Il69l4l77l Il82l Il98j | . 
Different and independent data of many observations, such as Wilkinson Microwave Anisotropy Probe, SNe la and 
BAO, were specially analyzed in order to i nvest igate the limitations on the intensity and form of the interaction in the 
dark sector. Some researchers also suggest |l92j |. that dynamical equilibrium of collapsing structures, such as clusters, 
will essentially depend on form and sign of the interaction between dark matter and dark energy. 



the problem is that during all the Universe history the two densities decay by different laws, so it is necessary to impose very strict 
limitations on their values in early Universe in order to make them be of comparable order nowadays 



28 



A. Model of Universe with time dependent cosmological constant 

The simplest example of the model with interacting DM and DE is the cosmological model with decaying vacuum. 
Actually the A(t)CDM cosmology represents one of the cases where parameter w for dark energy equals to —1. 

This model is based on the assumption that the dark energy is nothing that physical vacuum, and energy density 
of the latter is calculated o n th e curved space background with subtraction of renormalized energy density of physical 
vacuum in the flat space jl93j |. The resulting effective energy density of physical vacuum depends on space-time 
curvature and decays from high initial values in early Universe (at strong curvature) to almost zero magnitude in the 
present time. 

Due to the Bianchi identity, the decay of vacuum must be accompanied by creation or mass increase of dark energy 
particles, which is common property of the decaying vacuum models, or in more general case, for the models with 
interacting dark matter and dark energy. 

Now let us give the formulation of the A(t)CDM. In the present section we will use the system of units where the 

reduced Planck mass equals to unity: AIpi — {8nG) ' =1. 

For the case of flat Universe described by FRW-metric, the first Friedmann and the conservation equation can 
be presented in the form 

ptot = ptot + 3H{ptot + Plot) = 0, 

where ptot = p„i + p^i ^nd taking into account that p.,n = and p\ = —p\, the conservation equation takes the form: 

p„i + 'SHp,n = -pA- (137) 

The right hand side of the latter contains an additional term which plays role of a source — it is the decaying 
cosmological constant. The problem with such approach is in fact that we end up with the same number of equations, 
but acquire additional unknown function A{t). The above descr ibed approach to define the physical vacuum density, 
though intuitively clear, faces certain principal difficulties jl68j |. That is why a phenomenological approach prevails 
in the literature. Below we present a simple but exactly solvable model. Consider the case when A depends on time 
as [Zl: 

A = aH. (138) 

It is interesting to note that with the choice a (tott is the energy scale of the QCD vacuum condensation) the 

relation (|138|) gives the value of A which is very close to the observed one. 
For the considered components the first Friedmann equation takes the form 

p-y + A = 3H^. (139) 

The equations (|137p and (|139p together with the conservation equation 

P-t = Wp^ = (7- 1)^7' 

and the decay law for the cosmological constant (|137p completely determine the scale factor time evolution. Combining 
those equation, one obtains the evolution equation in the following form 

2H + 3-/H^ - (j-fH = 0. 

Under condition H > the solution for the latter equation reads 

a{t) = C[exp {<Tjt/2) - 1]^, 

where C is one of the two integration constants (recall that the equation for the scale factor is of second order). The 
second integration constant is determined from the condition a{t = 0) = 0. Using such solution, one can find the 
densities of matter (radiation) and time-dependent cosmological constant 



1 



(140) 



(141) 
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Let us now analyze the history of Universe expansion in the scenario under consideration. In the radiation-dominated 
epoch 7 = ?ii-|-l = 4/3 and, therefore, 



a{t) = C 



e3 



1 



For small values of time {at 1) and we recover the well-known dependence a cx i^/^. 



The densities p-,, (|140p and A (|141|) in the radiation dominated epoch transform into 

P-r ^ Pr = 



a^C^ 1 a^C^ 1 



3 a4 



, a^C^ 1 

3 3 



3 a2' 



In the linit a — > (t — > 0) one has 



Pr = 



3 a4 - 4^2 ; 



A 



(T^C^ 1 a 



3 a2 2i' 

Now consider the matter dominated era. In that case 7 = w + 1 = 1 and 



a(t) = C 



1 



For crt <C 1 we reproduce the standard law for matter evolution 

For the densities p^ and A in the matter dominated era one finds 
a^C^ 1 a2C3/2 1 



P7 ^' Pti 



A 



3 a3 



^2^3/2 1 



7,3/2 ^ 



3 3 a3/2- 

Note that in the limit of long times {at ^ 1) the scale factor growth exponentially as 
a{t) = Cet*. 

In the same limit the matter density tends to zero, and the time-dependent function A{t) transforms into the "real" 
cosmological constant. 

In the matter dominated epoch the Friedmann equation (jl39p can be presented in the form 
H{z) = Ho[l- + ^Ml + zf^ 
This expression can be used for calculation of the deceleration parameter (g — — l) 

ln^o{i + zf^ 



q{z) 



(142) 
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One thus finds for the current value of the deceleration parameter 
3 

q{z = 0) = -ilmo - 1. 

Therefore the accelerated expansion of Universe occurs under the following condition 
2 

^rnO < 

This condition is in fact satisfied for the observed value flmo ~ 0.23. From (|142p it follows that the transition from 
the decelerated expansion to the accelerated one took place at 

r 1 O 1^/3 

This value {z* sa 1.2) exceeds (being however of the same order O (1)) the corresponding value for SCM (z* 0.75) , 
which is the result of the matter production in the vacuum decay process. 

Therefore the model based on the time-dependent cosmological constant, where vacuum density linearly depends 
on Hubble parameter, appears to be quite competitive. It sufficiently accurately reproduces the "canonic" results, 
relative both to the radiation dominated and matter dominated eras. The present Universe expansion is accelerated 
according to the model under consideration. Verification of the model with the latest observational data obtained for 
SNla leads to the results (e.g., 0.27 < ^Imo < 0.37), that agree very well with the currently accepted estimates for 
the accelerated Universe expansion parameters. 



B. Interacting dark matter and dark energy 

One of the most interesting properties of dark matter is its possible interaction with dark energy. Although the 
most realistic models (SCM in particular) postulate that the dark matter and dark energy are uncoupled, there 
is no serious evidence to consider such assumption as a principle. In the present time many active researches are 
provided to investigate the possible consequences of such interaction [163, 166, 167, 169 177, 182, 192, 198]. As we 
have mentioned above, the interaction of dark components can at least soften some sharp cosmological problems, such 
as the coincidence problem for example. The dark energy density is approximately three times as high as the dark 
matter one. Such coincidence could be explained if dark matter was somehow sensitive to dark energy. 

Remark that the possibility of interaction between the dark energy in form of scalar field and dark matter lies in 
the basis of warm inflation model. Unlike the cold inflation scenario, the former does not assume that the scalar fleld 
in the inflation period is isolated (uncoupled) from other flelds. That is why instead of overcooled Universe in the 
inflation period, in the Universe model under discussion a certain quantity of radiation is always conserved, which is 
sufflciently remarkable to manifest in post-inflation dynamics. 

Interaction of the scalar field quantum, responsible for the infiation start - the inflaton - with other fields imply 
that its master equation contain s the terms describing the process of energy dissipation from inflaton system to 
other particles. Barera and Fang [l87j | initially assumed that consistent description of the inflaton field with energy 
dissipation requires the master equation in form of the Langeven equation, where the fluctuation-dissipation relation 
is present, which uniquely links the field fluctuations and dissipated energy. Such equation lies in the the basis of 
description of the warm inflation process. 

Interaction between the components in the Universe must be introduced in such a way that preserves the covariance 
of the energy-momentum tensor r.|"^, = 0, therefore T ^'^ .„ = —T^'^.,, ^ 0, where u^, is the 4-velocity. The 

[tot) ]U ' DM DE ' ' ^ 

conservation equations in that case take the form: 

u^T = ~u^T = -Q. (143) 
DM ;/i ^ DE ;m ^ ^ ' 

For FRW-metric the equations (|143p transform to: 

Pm + 3Hp,n = Q, (144) 

Pde+3Hp^J1 + w^^)^ -Q, (145) 

where pm and p^^^ are densities of dark matter and dark energy respectively, is the parameter of state equation 
for for dark energy, H = a/a is the Hubble parameter. 

„ r > ^ . f DE DM; 

^ < ^ ''''^'■Sy DM ^ DE 
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If Q < 0, then dark matter continuously decays into dark energy, if Q > then vice versa. Note that the equations 
(|144p and (|145p obey the conservation equation: 

Plot + iHptotiptot + Plot) = 0, (146) 

where ptot = Pde + Pm is the total energy density. 

The interaction between the dark matter and dark energy is effectively equivalent to modification of the state 
equation for the interacting components. Indeed, the equations (|144p can be transformed to the standard form of the 
conservation law written for uncoupled components: 

Pde + SiJPoB (1 + ^DE^.f! ) = 
Pi^A/ + = 0, 



where 

Q 



(147) 



play role of effective state equations for dark energy and dark matter respectively. 

As we do not kn ow t he nature of both the dark energy and dark matter, we cannot derive the coupling Q from 
the first principles '185"]. However it is clear from dimension analysis that this quantity must depend on the energy 
density of one of the dark component, or a combination of both components with dimension of energy density, times 
the quantity of inverse time dimension. For the latter it is natural to take the Hubble parameter H. 

Different forms of Q the most often used in the literature, are the following: 

Q^SH^Pm, Q = 3Hjp„,, Q = SHj{p^^ + p^J. (148) 

For example we consider the simplest case for such model 

Pm + iHpm = SHpm: 

Pde +3H{p^^ +PnE) = SHp,,,, 

where d is the dimensionless coupling constant. 
Integration of the latter equations yields 

p a-[3(i+«'i.B)+<5]. 

roE l-'nEo"' ' 

„ _ -^PpEO „-[3(l+tu„„)+i5l , f ^ I ^PnEO 



(149) 



Inserting the obtained densities into the first Friedmann equation (jl3p and transforming from scale factor to redshift, 
one gets for H'^{z) the following 



{l + zfHl r 



3(3t«BE +S) I i 
Note that in the considered case the deceleration parameter q{z) can be easily determined by the following 



, , 1+zdH^ ^ 



and direct substitution of (|150p results in explicit dependence of deceleration parameter on the redshift: 

, , 3 w^,nDE{Hi + w,,) + 6){i + z)^^^nE+^) + n,n + s/w^, 

qiz) — —I -\ -TT . (151 

2 3w^,nDEil + z)^^^oE+S) +n„, + S/w^, 

With Kjj^^ = — 1 and 5 = 0, this expression coincides with the deceleration parameter value obtained in SCM. Note 
that in the considered model appropriate choice of the coupling parameter value S can essentially soften or even get 
rid of the compatibility problem for the densities of dark matter and dark energy. Indeed, using the densities ()149p , 
one finds the following relation R = Pm/ Poe • 

R = : r +lRn + '- -] a'^'^oE+\ 



3w^E +S \ 3u;^_g 
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where Rq = Pmo/PoEo ratio of the densities in the present time. In SCM it is known to be i? '-^ a~^, which 

differs from the considered model on the S in exponent. 

Let us now consider the inverse problem, and instead of interaction we postulate the ratio 

— = f{a), (152) 

PnE 

where /(a) is arbitrary differentiable function of scale factor. Thus one obtains: 

Pr,i = PoEfia), (153) 

Pm 



roE r/ \ ' 

/(a) 

Pra= PoeI + PoEf'O; (155) 
Pm PmO-f 



(154) 

(155) 
(156) 



.OE f p ■ 

Inserting the expressions (|155p and (|153p into the equation (I144p we obtain: 

PdeI + PDEf'a. + 3Hp^^f ^ Q, (157) 

where prime denotes the differentiation with respect to the scale factor. Using here the expression for p^j^ , obtained 
from (|145p . one gets: 



1 + / V / 

The first Fricdmann equation takes on the form: 

= p^^ + p,n = Per, (159) 

where pcr is the critical density. Accordingly one can write down that 

f^.. - — = T^; (160) 

Per 1 + / 

^m = ^ = -^T- (161) 

Per 1 + 7 

Now, substituting (|16ip into (|158p . we finally obtain the expression for Q: 

Q = Hp^,n„, (l^ - iw^,^ = Hp„,n^, - iw,,^ . (162) 

It is worth noting that in the case when /(a) oc we obtain the expression 

g = Hp^n^^ (C - 3w^^) ^ Hp^^n^ {£_ - Sw^^) . (163) 

From the equation (I146P it follows that: 

d\n Ptot 
din a 

where 



-3{l+Weff), (164) 



Weff = = = —7 = il^^W^E, (165) 

Ptot PoE + Pm 1 + / 



and therefore 
rfln Ptot 



din a 

From the latter equation we obtain: 



'S{l + n^,w^,). (166) 



Ca ^ exp ^— y 3r2ogWj5£, dlna 



Ptot = Ca exp - / 3r2oEW^E c^lna (167) 
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where C is the integration constant determined from the requirement that ptot{a = 1) = ptot,o = 3Hq /{8ttG). Using 
the expression for ptoti one can easily obtain from the Friedmann equations the relation for E = H/Hq, which is used 
to verify cosmological models and establish restrictions on cosmological parameters. Expression for = ^l^^ptot 
and pm = ^mPtot can also easily obtained from the equations (|160p and (|167p . 
At last it is worth noting that one can obtain from the equation (|152p . 

/o = f{a = 1) = = _ = ^ 1. (168) 

Let it be now 

/(a) = /oa«, (169) 

where ^ is a constant and fo is defined above. Substituting the equation (|169|) into the equation (|167|) and requiring 
that ptotia = 1) = ptot,o, we can determine the integration constant and find that 

Ptot = Ptot,oa~^ [ilmO + - ^rno)a^] (170) 

Inserting this expression into the Friedmann equation, one finds: 



= = [f^^o + (1 - f^™o) a^r'^''^^ 

Hi (171) 

(1 + zf [ r!„,o + (1 - 0,„o) (1 + z)-^] ' ^ . 



Using the above obtained formulae, it is easy to find the following expression for the deceleration parameter in the 
considered model 

3 

9 = 1 + 2 ("'be^be + W'mf^m) 7 (172) 

where relative densities rij-,^ , the following form 

Q ^DEo Q _ (1 ~ ^deo)q^ (^^X\ 

"-"17,,„ + (1-I7„,ja?' ---17,,„ + (l-f],,Ja«' ^ ^ 

and we finally obtain 

^13 w^^n^^„(l + z)^ 

2^21-f],,„+17,,„(l + z)«' ^ ^ 

As the dark energy density is not constant in the considered model, unlike the SCM, it facilitates solution of the 
coincidence problem under condition ^ < 3. In the model under consideration the accelerated expansion phase starts 
at the redshift z* , which is defined by the relation q{z*) = 0, and it equals to 



1 ^r>EO 



(1 + 3Wbe)^^ceo 



1. (175) 



This relation naturally generalizes the expression ((58)) . obtained for SCM. The difference of the value (|175p from the 
corresponding one for SCM depends on magnitude of the distinction of the parameters ^ and from 3 and -1 
respectively. In the case of zero difference the quantities (|55|) and (11751) coincide. Note that at z — )■ oo Universe 
expanded with deceleration g — > i, the case z — ^ — 1 responds to g — )■ i — l^^^. Therefore, as it was expected, 
dynamics of the considered model is asymptotically (at z — > oo and z — )■ — 1) identical to the case of two-component 
Universe, filled by non-relativistic matter and dark energy with the state equation p — w^^^p. 



C. Cosmological models with the new type of interaction 

In the present subsection we consider one more type of interaction Q, [t^ which change its sign, i.e. direction of 
energy transfer, when the expansion changes from decelerated regime to accelerated one and vice versa. 
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Some recently appeared papers |79| make attempt to determine both the very possibihty of interaction in the dark 
sector and specify its form and sign, basing solely on analysis o observational data. The analysis splits all available 
set of redshift data z into some parts, within which the function 8(z) — Q/{ZH) is considered to be constant. This 
analysis allowed to establish that the function 5{z) the most likely takes zero values {5 = 0) in the interval of redshifts 
0.45 z < 0.9. It turns out that this remarkable result produces new problems. Indeed, as we have already mentioned, 
the researchers mostly consider the interaction of the form ()148|) : Q = 3Hjpm, Q = SHjp,-,^^, Q ~ 'dH^{pm + Poe) 
and for the model under consideration it is always either positive or negative defined, and therefore it cannot change its 
sign. Sign changes are possible only in the case when 7 = /(t), which can change the sign of Q or Q = iH{apm+PPoE)^ 
where a and j3 have opposite signs. 

As authors of [79{ mention, solution of such problem requires to introduce a new type of interaction, which can 
change its sign during the evolution of Universe. 

An interaction Q of such type was proposed in the paper [t^ , where its cosmological consequences were considered. 
As it was noted in the above cited paper, the redshift interval, where the function 5{z) must change its sign, includes the 
transition point when the Universe expansion stopped to decelerate and started the acceleration (see ([58])). Therefore 
the simplest type of the interaction which can explain the above mentioned property is the case when the source Q 
is proportional to the deceleration parameter q : 

Q = q{ap + 3l3Hp) (176) 

where a and /3 are dimensionless constants, and the sign of Q will change with the transition of universe from the 
decelerated expansion stage {q > 0) to the accelerated one {q <0). The authors also consider the case 

Q^q{apm+3l3Hp,n), (177) 
Q^qiaptot+i^Hptot), (178) 
Q = qiap^,+3l3Hp^,). (179) 

The paper [zl] considers a model of Universe with decaying cosmological constant 

PA = -Q ■ 

The Friedmann and Raychaudhuri equation take thus the form 
2 ^2 

= yPtot = y(pA+P,n), (180) 

H = -{Ptot + Plot) ^ jPm, (181) 

where = SnG. Following the paper ^], in the succeeding subsections we consider cosmological models with 
interaction of the type ([T771) - p7g)) . 



1. Case Q = q{oipm + iPHpm.) 

For the beginning we consider the case when the interaction takes the form (I177P and insert such expression into 
the conservation equation (I144p , resulting in the following 

= • 3ilp„ . (182) 

1 — aq 

Substituting the obtained expression into the equation (|177p . we finally get 

Q = ■ SqHp,-,. (183) 

1 — aq 

From the equation (|18ip one obtains 

Pm^-^H- (184) 
Inserting it into the equation (|182p one finds that 

ii = -SHH , (185) 

1 — aq 
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Thus we obtained the second order differential equation for the function H(t). Now transform from the time derivative 
to differentiation with respect to the scale factor (denoted by the prime '), then the equation (|185p takes on the form 



aH'- 



H 



(186) 



1 — aq 

This expression represents a second order differential expression for the function H{a). Note that the deceleration 
parameter 

is also function of H and H' and, except the case a 7^ 0, the equation has no exact solution and it represents 
a transcendental differential equation of second order. That is why we consider solely the case a — 0. Thus the 
interaction (|177p takes the form 

Q = 3l3qHp,n- 

With a = Q the solution (|186p can be presented in the form 
H{a) = C12 [3Cii(l + /3) - (2 + 3^) a-'^'^^+^^ ^ 1/(2+3/3) 



(187) 



where Cn and C'12 are the integration constants determined below. We find the relative density of dark matter as 
the following 



2H 



2aH' 



3iJ2 3i72 3ij • 
Inserting the equation (|187p into (|188p . one gets 
2(l + /3) 

" 2 + 3;3-3Cii(l + /3) a3(i+« • 

With the requirements fim(a = 1) = r2„io and H{a = 1) = Hq, the integration constants take the form 
^ 0„o(2 + 3/3)-2(l + /?) 

C12 = i/o[3Cii(l + /3)-(2 + 3/3)]-^/('+=^^^ 
Substitution of the expressions (|190p and (|191|) into the equation ()187|) finally gives the result 



E = 



H 
Ho 



3/3 



2(1 + /3) 



l-(l + z) 



3(l+/3) 



1/(2+3/3) 



(188) 

(189) 

(190) 
(191) 

(192) 



The model contains two free parameters flmo and f3. We note that if /? = 0, then the equation (|192p reduces to 

1/2 

i?(z) — [n„io{l + z)^ + (1 — f^mo)] , which is equivalent to ACDM model. Using the relation 
(1 + z) d f 1 \ 

one finds the dependence of deceleration parameter on the redshift in the considered model 



3 (l + z)^(i+^) 

'ZW = -1 + ^^mO ^(2+3/3) ■ 

The effective parameter of the state equation is known to equal 

Ptot (2g - 1) 
WeS = = 5 ■ 

Ptot o 



(193) 



The figure m presents the plots for dependence of some cosmological parameters on the redshift z. The free parameters 
and /3 were chosen to provide the best agreement with observations. One can find that in the considered model 
the transition from the decelerated expansion (q > 0) to the accelerated one {q < 0) took place at Zt = 0.7489, the 
parameter /3 is negative and therefore dark matter decays into dark energy when z > zt, and vice versa at z < zj. 
The Universe is interaction- free in dark sector at Zf. 
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FIG. 4: Qm, ^A, q and WcS as function of the redshift z at Q.mo = 0.2738 and /3 = —0.010 in tiie case Q — SfiqHp, 



2. Case Q = q{aptot + ZPHptot) 

Now we consider the case (|178p . and proceeding completely analogous the above considered case one obtains 
3gi/3 



Q 



2a-^ + 3/3 



(194) 



Inserting the equations (|184p and (|194p into (|144l) . and transforming to differentiation with respect to the scale factor, 
one finds 



aH" + + (4 + iaq) H' + ^ = . 

H Za 



(195) 



As in the previous case, we obtained again the differential equation of second order for the function H{a). Exact 
solution exists only in the case a — : 



H{a) = C22 • a-3(2-30+'-i)/8 . (^^3r,/2 ^ ^^^j 



1/2 



where C21, C22 are integration constants and ri = ■\/4 + /3 (4 + 9/3). Inserting (|196p into (|188p . we get 



2-3/3 + 



2C21 



+ C21 



1 n 



The integration constants are determined as usual from the condition r2,„(a = 1) = ClmOi H[a = 1) = Hq : 



C21 = -1 
We finally get 
H 



2ri 



2 - 3/3 - 4f7,„o + n 



, C22 — (1 + C21) 



^1/2 



E= — ^ (1 + ^)3{2-3/3+ri)/8 



(1 + z)-3-i/2 + C21 



I + C21 



1/2 



(196) 



(197) 



(198) 



(199) 



In the considered model there are also two free parameters f2mO and /3. Using the condition < < 1 with a — 0, 
from the equation (|197l) one gets /3 > 0. The best agreement of the model under consideration with observational data 
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occurs with flmo = 0.2701 and /3 = 0.0. It means that the considered interaction model gives worse agreement with 
observations than AC DM. More detailed discussion an interested reader finds in the paper [t^ by the author of the 
considered model. The transition form the decelerated expansion phase {q > 0) to the accelerated one {q < 0) occurs 
at zt = 0.7549. 




0.75 




6 7 



FIG. 5: The same as on Fig. |4l but in the case of interaction of the form Q — SPqHptot under condition /3 > 0. 



3. The case of Q = q{apA + SPHpA 



For the conclusion we consider the case (|179p . Following the same procedure as in the two preceding case, one 
obtains 



Q 



3/3gi?PA 



I + aq 

With the equation (1184^ one has 

PA = - = ^ (SH^ + 2H 

Therefore the equation for the Hubble parameter in terms of the scale factor takes the form: 



aH" 



12 



3Pq 



H' 



9l3qH 
2a(l + aq) 



= 0. 



H \ 1 + aq^ 

Exact solution can be obtained in the case, Q — 3f3qHp\, namely 

H{a) ^ C32 • a-3(2-5/3+'-.)/[4(2-3/5)] . (^^3r,/2 ^ ^31)'^^'"'"^ , 



(200) 



(201) 



(202) 



(203) 



where C31, C32 are the integration constants, and r2 = a/(2 — /3) — \ 2 ~ (3 \ . Inserting (|203p into (|188|) . we get 



1 

2 (2 - 3^) 



2-5/3 + 



2C31 



-1^2 



(204) 
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FIG. 6: The same as Fig. [H but for the case Q — SPqHpA- 



Assuming that flm{a = 1) = ^mo and H{a = 1) = Ho, wc can write 

2 '"2 ^ rr . ^ ^l/(3/3-2) 



C31 = -1 

and finally get 
H 



2 - 5/3 + r2 + 2n,nO (3/3 - 2) ' 



C32 — Ho{l + C31) 



E 



Ho 



(1 + 2)3(2-5/J+r2)/[4(2-3/J)] 



(1 + + C31 



I + C31 



1/(2-3/3) 



(205) 



(206) 



As before, the model has two free parameters flma, P- The maximum plausibility method for the free parameters of the 
considered model gives the result [T^I ilmo = 0.2717, /3 = 0.0136. Unlike the two preceding models, the observational 
data analyzed in [78|, give evidence in favor of /3 > 0. More detailed discussion the interested reader can find in the 
paper (tI] of the author of the model. 



The plot [B] displays the dependence for deceleration parameter and effective parameter of the state equation = 
Plot I Plot — (2(7— 1)/3 as functions of the redshift z with the parameters obtained by the maximum plausibility method. 
It is easy to show that the transition from the decelerated expansion phase (g > 0) to the accelerated one (g < 0) 
takes place at zt — 0.7398. As the parameter /3 obtained from observations satisfies /3 > 0, then dark energy decays 
into dark matter [Q > 0) for z > zt, and vice versa (Q < 0) for z < zt- 



D. Statefinders for the interacting dark energy 

In this part of the review, we discuss recently introduced "statefinder parameters" [o^, that inclu de the third 
derivative of the cosmic scale factor, are useful tools to characterize interacting quitessence models [l28j |. 

Accelerated expansion of Universe is normally described by the deceleration parameter q = —d/{aH'^). As was 
mentioned above, current value of the Universe deceleration parameter q is negative, however its amplitude is difficult 
to determine from the observations. Since the accelerated dynamics of Universe is predicted by many models, they 
can be used to obtain an additional information on q. 

In particular, an example of such models is given by the cosmological models where the dominant components - 
dark matter and dark energy - interact with each other. The models, where the dominating components do not evolve 
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independently and interact with each other instead, are of special interest because, as was mentioned above, they can 
facilitate or even completely solve the "coincidence problem" . 

Let us consider the Universe filled by two components of non-relativistic matter (indexed by to) with negligibly 
small pressure p„i ^ Pm and dark energy (indexed by x) with the state equation = wpx, where w < 0. The dark 
energy decays into dark matter according to 



Pm 
Px 



3i?Pm — Q , 
3H(l + w)p, 



Q, 



where Q > measures the strength of the interaction. For further convenience we will write it as Q 
n is a new variable with dimension of pressure. 

The Einstein equations for spatially fiat Universe have the form 



H 



\ttG 



-P , 



{P+Px) : 



(207) 
-3niJ, where 

(208) 
(209) 



where p = Pm + Px is the total energy density. The quantity H is connected by the deceleration parameter q by the 
relation q — —1 — [H / H"^) = (1 + 3it;51a;)/2, where 0,^ = ^ ttG px / [iH"^) ■ Evidently the deceleration parameter does 
not depend on interaction between the components. Nevertheless, differentiation of H yields 



H 



9 



Px 

p 



, ^ Pa; n W Px 
W(l+W) W XTT — 

p p SH p 



(210) 



In contrast with H and H , the second order derivative H in fact depends on the interaction of the components. 
Therefore, in order to distinguish models with different types of interaction, or between interacting and free models, 
we need the cosmology dynamical description involving the parameters explicitly dependent on H . In the papers 
and [i^ Sahni and Alam introduced the pair of parameters (the so-called " statefinders" ) , which seem to be promising 
candidates for investigation of distinction between the models with interacting components. 



1 



3(9-^) 



In the considered context of interacting components they take on the form 



9 w 



2 1 



w - 



n 

Px 



3wH 



where k = Pm/ Px, and 

n 

s — 1 + w 

Px 



iwH 



(211) 



(212) 



(213) 



For the interaction- free models with 11 = 0, the statefinders reduce to the expressions that were studied in [95| and [96|. 
Let us show how the statefinders can be useful for analysis of cosmological models, where the dominant components 
interact with each other. Note that the third order derivative of the scale factor is necessary for description of any 
variations in the general state equation for the cosmic medium [l9| . It is obvious from the general relation j96| 



1 



9p + P P 



2 P p 

where P is the total density of cosmic medium which in our case leads to the relation P px- Since 



P 
P 



P 
P 



(214) 



(215) 



it is evident that the interaction term in the relation P px — wpx 
affect the time dependence of the common state parameter P/p. 



WPa 



in accordance with (|207p will additionally 
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Models with interacting components suggest a dynamical approach to the coincidence problem. The principal 
quantity in that context is the density ratio k, which was defined above in expression l|212p . Strict solution of the 
coincidence problem requires that the density ratio should be of order of unity on late stages of Universe evolution. 
At least it has to vary sufficiently slowly for time interval of order H~^. The ratio k is governed by the equation (|207p 



-3H 



Pm Px 



(216) 



Now let us consider the statcfinders for different solutions of that equation. 



E. Scaling solutions 



In the paper [98| it was shown by the authors that the so-called scaling solutions, i.e. the solutions of the form 
Pm/px oc a~^, where ^ is constant parameter, lying in the range [0,3], can be calculated in the case when when 
dark energy decays into dark matter (see eq. ()207|) ). The solutions of such type make interest because they facilitate 
solution of the coincidence problem [99[. Indeed, the model with ^ = 3 reduces to ACDM one for zz; = — 1 and H = . 
For ^ = the Universe dynamics becomes stable with k = const, and the coincidence problem does not appear [l00j |. 
Thus for the parameter values ^ < 3 the coincidence problem weakens. In this scheme with w = const, one can easily 
see that the interaction which generates the scaling solutions, can be calculated as the following 

n / n «o(i + ^)^ 

where kq denotes the current ratio of the energy densities and z — [a^/a) — 1 is the redshift. Substituting this 
expression into the equations (|212p and (|213p . one finds the statefinder parameters: 



9 

r = 1 + - 



2 1 + Ko(l + z)« 



w — \ w -\ 



Ko(l + z)^ 

1 + Ko(l + 



and 



S = 1 + Ui 



«q(1 + ^)^ 

1 + K0(1 + Z)«' 



(218) 



(219) 



Figure [7] presents the dependence r(s) for different values of ^. The less is the value ^, the lower is the corresponding 
curve on s-r diagram and the less actual is the coincidence problem. Remark that those curves are qualitatively 
identical to those drawn for the models with interaction- free components (see FiglS] 95] ) . Compare also that ACDM 
model (n = 0, u> = — 1) corresponds to the point s = 0, r = 1 (not shown on the figure). The present values rp 
and So of the statefinder parameters are important both from the observational point of view and for the purpose of 
distinction between different models. For the scaling models one has 



, 9 
'^0 = 1 + :^- 



2 1 -t- Ko 
Having in mind that 

1 1 -I- Ko + 3w 
*=2 1 + Ko 
and introducing 

9oA = go (w = -1) = 



and 



So = 1 



Ko 



1 + Ko 



12 



_Ko^3 



Kq 



= 1 



2 1 + Ko 2 1 -I- Ko 

we can classify different models by its dependence so(go 



90A , 

), which is 



<7o - goA + 



(1 + goA) 



(220) 



(221) 



(222) 



(223) 



The first part of the bracket in the righthand side describes deviation from w = —1, and the second responds for 
deviation from ACDM value ^ = 3. For models with w = —\, which have the same deceleration parameter go = goAi 

we obtain sq = f (| ~ + ^oa) • Of course, ^ = 3 corresponds to ACDM model with sq = 0. If the conditions 

Ko = 3/7 and ^ = 1 hold, then so (C = 1) = -0.2, while for ^ = one has Sq (^ O) = -0.3. 

Analogous considerations are valid also for other values of w. Therefore the parameter so enable us to distinguish 
between different scaling models, which have the same deceleration parameter. 
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FIG. 7: The curves r(s) are drawn for the redshift range [0, 6] (from left to right) with w = —0.95 and kq — 3/7 for 3 different 
values of the § parameter: (a) 2.5; (b) 1.5; (c) 0.5. 
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Statefinder s 

FIG. 8: The Statefinder pair (r, s) is shown for different forms of dark energy. In quintessence (Q) models {w =constant 7^ — 1) 
the value of s remains fixed at s = 1 + u; while the value of r asymptotically declines to r{t S> to) — 1 + ^(1 + w). Two 
models of quintessence corresponding to wq — —0.25, —0.5 are shown. Quintessence (K) models are presented by a scalar 
field (quintessence) rolling down the potential V{<l)) oc tf)'" with a = 2,4. ACDM (r = l,s = 0) and SCDM in the absence 
of Z (r = 1, s = 1) are the fixed points of the system. The hatched region is disallowed in quintessence models and in the 
quintessence model which we consider. The filled circles show the current values of the statefinder pair (r, s) for the Q and K 
models with flmo = 0.3. 



F. Concluding remarks 

The statefinder parameters introduced in [95| and [11), as was assumed, became useful tools for verification of 
models with interacting components, which in turn solve or weaken the coincidence problem — a stone of obstacle 
for many models with late acceleration. While the deceleration parameter does not depend on interaction between 
the dark energy and dark matter, the state finders (r, s) contain the dependence explicitly. 
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IX. HOLOGRAPHIC DYNAMICS: ENTROPIC ACCELERATION 



In the previous section we presented the description of Universe dynamics based on general relativity with some 
generalizations. The present section is intended to show that there exists a principally independent approach, which 
enables us both to reproduce all achievements of the traditional description and resolve a number of problems the 
latter faced to. 

The traditional point of view assumed that space filling fields constitute the dominant part of degrees of freedom in 
our World. However it became clear soon that such estimate much hardened the development of the quantum gravity 
theory: for the latter to make sense one had to cutoff on small distances all the integrals appeared in the theory. 
Consequently our World was described on a three-dimensional discrete lattice with cell size of order of Planck length. 

Recently some physicists came up with even more radical point of view: complete description of Nature required only 
two-dimensional lattice, situated on space edge of our \yorld instead o f the three-dimensional one. Such approach is 
based on so- called " holographic principle" ^M, [iM [HI ^M^M ■ 

The term comes from the optical holography, 
which represent nothing but two-dimensional recording of three-dimensional objects. The holographic principle is 
composed of the two main statements: 

1. all information contained in some region of space can be "recorded" (presented) on boundary of that region; 

2. the theory contains at most one degree of freedom per Planck area on boundaries of the considered space region 

N < (224) 

Therefore central place in the holographic principle is occupied by the assumption that all information about the 
Universe can be coded on some two-dimensional surface — the holographic screen. Such approach leads to possibly 
new interpretation of cosmological acceleration and to completely novel concept of gravity. The gravitation is now 
defined as entropic force generated by variation of information connected to positions of material bodies. 

Let us consider a small piece of holographic screen and approaching particle of mass m. According to the holographic 
principle, the particle affects the information amount (and thus on entropy), which is stored on the screen. It is 
naturally to assume that the entropy increment near the screen is linear on the displacement Ax, 

777 n 

AS^2TrkB — Ax. (225) 

h 

The factor 2tt is introduced for convenience reasons, to be clear below. To make clear why that quantity is proportional 
to mass, let us imagine that the particle splits on two or more particles of lower mass. Each of them then carries 
its own entropy increment with displacement on Aa;. As both entropy and mass are additive function, then they 
are proportional to each other. According to the first principle of thermodynamics, the entropic force, caused by 
information variation, reads 

FAx = TAS. (226) 

If the entropy gradient is given, or calculated using the relation (I225p . and the screen temperature is fixed, then one 
can find the entropy force. As is known, an observer, moving with acceleration a, effectively feels the so-called Unruh 
temperature jl44j 

in , ^ 

keT = a. 227 

In c 

Let us assume that the total energy of the system equals E. A simple assumption follows, that it is uniformly 
distributed on all N bits of information on the holographic screen. The temperature is then defined as average energy 
per bit 

E = ^NksT. (228) 

The formulae (|224p - (|228p enable us to formulate the holographic dynamics, and in particular dynamics of Universe, 
without any notion of gravity. 

The above presented ideology is essentially based on successful description of black holes physics [I5ll4l54 |. On 



the first face, there is nothing common between the extremely diluted Universe with density p ^ 10 g/cm? and 
"typical" black hole of stellar mass with critical density p ~ IQ'^'^ g / cm^ [M = IQMq). But the situation sharply 
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changes if we transit to more massive black holes. As we have seen above, the black hole radius scales cx M, 
therefore for the black hole critical density one has p oc M~^. Super- massive black holes situated in the galactic nuclei 
reach masses of order lO^^M©, and their gravitational radius (r^ ~ 3 x IQ^^cm) five times exceeds the dimensions 
of Solar system. The critical density in the latter case equals p ^ IQ^^ g / cm? , which is one order of magnitude less 
than the room air density (1.3 x g / cw}) . Let us consider, how close are the parameters of our observed Universe 
to those of black hole. We estimate mass of the observed Universe taking the Hubble radius for its radial dimension 
The mass contained inside the Hubble sphere thus reads 

47r 

M^mv = y^hP- (229) 

In the case of fiat space , the density p of the observable part of Universe enclosed by Hubble sphere can be reasonably 
replaced by the critical one p — 3H'^ / (SttG), and therefore the formula (|229p takes the form 

M ■ 

Inserting the obtained mass value into the expression for gravitational radius of the Universe, Vg = 2GMuniv, one 
finds 

rg = RH- (230) 

The obtained result for many reasons is nothing but a rough estimate, however it obviously favors the above cited 
black hole argumentation in application to the observable Universe. Remark that in the case of the Sun the ratio of 
the physical radius Rq = 695 500fcm to the gravitational one rgQ = 3km counts more than five orders of magnitude. 

Now let us identify the holographic screen with the Hubble sphere of radius R = (which is valid for fiat 

Universe) and we try, making use of the holographic principle, to reproduce the Fricdmann equations, using neither 
Einstein equations nor Newtonian mechanics. The holographic screen spans the area A ~ AttR^ and carries the 
(maximum) information N = 47r i?^/Lpj bits. The change of information quantity dN for the time period di, caused 
by the Universe expansion i? — >■ i? + dR, equals 

dA 8ttR 
dN=-^^ -j^dR. 

^pi ^Pi 

Here we use the system of units c — — I- Variation of Hubble radius leads to change in Hawking temperature 

From the equipartition law it follows that 

dE = ]-NdT + ]-TdN = -^dR = ^, (231) 
Z Z PI 

[L^i = ■^) ■ The quantity dR can be presented in the form 

dR = -HHR^dt. (232) 

From the other hand, the energy fiow through the Hubble sphere can be calculated if the energy-momentum tensor 
is given for the substance filling the Universe. Considering the latter substance as an ideal liquid, (see chapter 2) and 
using the relation T^^ = {p -\- p)ui^iUi, +pg^v: we obtain 

dE A{p + p)dt. (233) 

Equating to ((255)) accounting (gM]), one finds 

H = -47rG {p + p). 



as we have repeatedly mentioned above, deviations from the flatness are very small 
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It is well known that the system of equations 

H = -4TrG{p + p), p + 3H{p + p)= 0. (234) 



is equivalent to the system (|T3|)-([T4t. Thus the declared aim is achieved. We remark that already in the year 1995 a 
more general problem was solved [119|: the Einstein equations were obtained from the thermodynamic considerations. 
This important result followed from the fact that entropy was proportional to horizon area and from the assumption 
that the relation (|227|) holds for any accelerated moving observer situated inside her own casual horizon, and T stands 
for Unruh temperature. 

Derivation of the Friedmann equations from the holographic principle is of course an important result, but by 
itself it represents nothing but reproduction of something already well-known. A natural question is whether one can 
develop a novel approach to description of Universe dynamic basing on the holographic principle? If yes, than is it 
possible to overcome the unsolvable difhculties of the traditional approach in frames of the new one? Let us start 
from the logic scheme which the holographic dynamics of Universe is based on. 

The standard procedure to develop any field theory involves the formulation of equations of motion obtained from 
some effective action corresponding to the energy scale characteristic to the phenomenon to be described by the theory. 
General relativity theory is not an exception from the rule: equations of motion in general relativity are obtained by 
variation of effective action with respect to dynamical variables present in it. The standard approach to evaluate the 
integrals that arise in variations of the action is the integration by parts, which produces the terms corresponding 
to boundaries of the region under consideration. According to such approach a quite reasonable assumption is 
made that on boundary of the considered region, far from the field sources, all the dynamical variables, be it either 
fields, or gravitational potentials or connectivity components (as in general relativity), turn to zero. In holographic 
physics instead the surface is the place where the main action takes goes on. The natural conclusion follows that the 
contribution of surface terms now should be necessarily accounted. 

Let us show that consideration of the boundary term in the Einstein-Hilbert action is equivalent to inclusion of 
non-zero energy-momentum tensor in the standard Einstein equations. The action for gravitational field takes the 
form ^ 



d'^xy^R. (235) 



167rG, 

Variation of that action with respect to metrics g^i, in a compact region fl, one obtains 



6 ( d^xy/^R 
Jn 



(236) 



where = q°^V^ V^. Full derivative in (|235]) can be presented as the contribution of the region boundary dfl, where 
those contribution do not cancel. In order to obtain the Einstein equations from the least action principle without 
boundary terms, the Einstein-Hilbert action should be complemented by a functional intended to compensate the 
contribution of full derivative in (|235p . We denote it as SBoundary[g]-- Then the full action takes on the form 

S = — [ d'^X^/^R-\' SBoundary[g\ + S Source, (237) 

iDTTLT 



n 



where we included possible sources Ssource of gravitational field, coupled with the material fields. The functional 

S Boundary [g], thus takcs the form 



SS 



Boundary 



[g] = - Jj'xV^^g^-'V'Sg^, - V^V^%,) , (238) 



In the context of holographic physics where the boundaries play crucial role, the functional Ssoundaryig] can be 
treated as action for holographic dark energy. The Einstein equations including the contributions of material fields 



t We use the signature {+, -, -, -), and definition R"^^^ = O^Fg^ - d^F^^^ + V.r^^ - Kx^^i^r, = Kb/' ^ = 9"" 
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and boundary terms have the fohowing form 

^ 2 3^"^^ ^ ^'^^ Source + Boundary) ' 

rpi-iii 2 SSsource rpiii> 2 SSBoundary /"OQQA 

-'-Source ~ i — - 77 ' ^Boundary ~ i — - 71 ' 

Remark that in the flat space case, when 17^,^ — 0, the boundary action SBoundary\g\ — 0, and Einstein equations take 
the standard form. We also note that such modification of Einstein equations preserves the general structure of the 
equations, which, as we will see below, guarantees the conservation of the Friedmann equation structure. 

Inclusion of the additional terms in the Einstein equations requires the so-called "holographic correction" of the 

Friedmann equations^ The easiest way to perform the correction is to insert the entropic force into the second 

Friedmann equation |l20l |. The entropic term structure can be guessed from the dimensional considerations 



a 47rG , , , , 

- = -— (p + 3p) + -^. (240) 

a, o 

Here ac—F^/m is the acceleration caused by the entropic force and Lb is the scale length, determined by the holographic 
screen position. If one takes the Hubble sphere as the holographic screen, then the scale length coincides with the 
Hubble radius: if, — Rh = H~^. The entropic acceleration magnitude can be expressed through the holographic 
screen temperature Tf, 

ae = ^=T~l = 2^rfc. (241) 
m l\x m 

Does the considered scheme agree with the main observational data at least? We have shown above in frames of 
SCM that absolute value of the cosmological acceleration of the Hubble sphere equals F ~ 4 x 10~^°msec~^. Let us 
estimate that magnitude in the holographic approach. 

It is naturally to assign the Unruh temperature to the holographic screen Tt, = Tjj- Let us link it with the Hawking 
radiation temperature Tu, which is equal to 

= STrfcsGA/ " 2^' ^^^^^ 
where g is the free fall acceleration on the Hubble sphere. Comparing the latter expression with Unruh temperature 

one gets 

TH^Tu{a = g). (244) 

Thus the Unruh temperature coincides with that for the Hawking radiation, but depends on the reference frame 
acceleration instead of surface gravitation (free fall acceleration) . According to the equivalence principle, the free fall 
acceleration on the Hubble surface is equivalent to acceleration of a reference frame and we can relate it to each other. 
Treating the Hubble sphere as an analogue of events horizon for a black hole, we can present the formula ()242p in the 
following form 

Th = AiL^Sx 10-30if. (245) 
Kb 27r 

Comparing (|245p and ([243]), one finds that the Hubble sphere acceleration an equals 

an = i? ~ IQ-Vsec-^ (246) 

The latter result is in convincing qualitative agreement with the one obtained in SCM. 

Certainly, quality of a novel model is not defined by its ability to reproduce the results established by its predecessors 
— it is nothing but minimum task. Its main purpose and virtue is to solve the problems untouchable by former models. 
The holographic dynamics is very promising from that point of view as it gives hope to solve both main challenges of 
the SCM — they are cosmological constant magnitude and coincidence problem. 
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Let us start from the former. The cosmological constant problem stands for the huge disagreement ^ between the 
observed density of dark energy in form of the cosmofogical constant and its "expected" value. The expectations are 
based on rather natural assumptions on the cutoff parameter for the integral, representing the zero fluctuations of 
vacuum. The holographic principle enables us to replace the "natural assumptions" by more thorough quantitative 
estimates. 



A. Universe with holographic dark energy 

In any effective quantum field theory defined in space region with typical length scale L and using the ultraviolet 
cutoff A, the system entropy takes the form S cx K^L^ . For instance, fermions, placed in nodes of space lattice with 
characteristic size L and period A~^, occupy one of the number 2'^^^^ states. Therefore entr opy of such system is 
S cx A'^L'^. According to the holographic principle, this quantity must obey the inequality |l2l| 

L^K^ < Sbh ^ \^ = ^LHdli, (247) 



where Sbh is the black hole entropy and Abh stands for its event horizon area, which in the simplest case coincides 
with the surface of sphere with radius L. This reasoning shows that magnitude of the infrared cutoff cannot be chosen 
independently of the ultraviolet one. 

Thus we formulate the important result |l2l| : in frames of holographic dynamics the infrared cutoff magnitude 
is strictly linked to the ultraviolet one. In other words, physical properties on small UV-scales depends on physical 
parameters on large IR-scales. In particular, if the inequality ()247|) holds, then one gets 

L^A-^Mli. (248) 

Effective field theories with UV-cutoff (j248p obviously include many states with the gravitational radius exceeding 
the region where the theory was initially defined. In other words, for arbitrary cutoff parameter one can find sufficiently 
large volume where entropy in an effective field theory exceeds the Bekenstein limit. In order check it let us remark 
that usually an effective quantum field theory must be able to describe the system under temperature T < A. While 
T ^ 1/L, the system has thermal energy M ^ L^T* and entropy M ~ L^T^ . The condition (|247p is satisfied at 

1/3 

T < i^Mpi/V) , which responds to gravitational radius Vg ^ L{LMpi) ^ L. In order to get rid of that difficulty an 
even more strict limitation [l2ll | is imposed on the infrared cutoff L ^ A^^, which excludes all the states localized 
within limits of their gravitational radii. Taking into account all above mentioned and using the expression (|249p 

Pvac^^, (249) 
the condition (|247p can be presented in the form 

L^PA < LMli = 2Mbh, (250) 

where Mbh is the black hole mass with gravitational radius L. Thus total energy contained in the region of dimensions 
L cannot exceed the mass of a black hole of roughly the same size. This result agrees surprisingly well with the 
expression (j230l) provided the IR-cutoff scale is identified with the Hubble radius H~^, besides that such choice is 
well motivated in cosmology context. The quantity pA is commonly called the holographic dark energy. 

In the cosmological aspect under interest the link between the small and large scales can be obtained from rather 
natural condition, that total energy contained within a volume of dimension L cannot exceed mass of a black hole of 
the same size: 

LVa < Mbh - iMpz- (251) 

Here p\ is the energy density. Would this inequality be violated, the Universe was composed of black holes solely. 
Applying this relation to the Universe as whole it is naturally to identify the IR-scale with the Hubble radius H~^. 
Then for the upper bound of the energy density one finds 

PA L-'^Mli - H'^Mli. (252) 



about 120 orders of magnitude !!! 
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The quantity pA is usually called the holographic dark energy. We will below denote its density as poE- Accounting 
that 

Mpi ~ 1.2 X lO^^GeF; Hq ~ 1.6 x IQ-^'^GeV, 
one finds 

PDB - IQ-^'^GeF''. (253) 

This expression will agrees with the observed value of the dark energy density ~ 3 x IQ-^'^GeV'^. Therefore the 
holographic dynamics is free from the cosmological constant problem. 

Though the obtained dark energy density corresponds to the correct one, the choice of the IR-cutofF scale equal to 
the Hubble radius a problem arises connected wi th th e state equation problem: the holographic dark energy cannot 
provide the accelerated expansion in that case (l58j . It can be easily seen from the following simple arguments. 
Consider the Universe composed of holographic dark energy with density given by the relation (I252p and of normal 
matter component. In that case p — pA + PM- From the first Friedmann equation it follows that p ex. . If pA oc , 
then dynamical behavior of holographic dark energy coincides with that of normal matter, thus the accelerated 
expansion is impossible. If the Universe was initially dominated by matter, then the dark energy defined by (j252p is 
in fact a tracker solution, as it reproduces dynamics of the dominating component. The presence of tracker solutions 
in frames of the holographic models gives a hope to solve the coincidence problem, but it contradicts to the observed 
dynamics of Universe. 

In order to produce the accelerated expansion of Universe in frames of holographic dark energy model we will try 
to use the IR-cutofF spatial scale different from the Hubble radius. First which comes in mind is to replace the Hubble 
radius by the particle horizon Rp: 



* dt r da 

— = a / 5 

a Jo Ha^ 



Rp^a I —^a I -j^. (254) 



Unfortunately such substitution cannot give the desired result. It can be made clear again using the initially matter- 
dominated Universe model. Due to the causality principle the gravity influence cannot extend to the regions situated 
on distances exceeding the causality horizon size. Therefore the vacuum energy density present in the Friedmann 
equations cannot be an arbitrary function of L. In the matter dominated era the causality (particle) horizon scales 
as a^/^, therefore p^j^{L) ~ a~^. As p{t) ~ a{t)~^'-^'^'^'\ then it immediately follows that for the case of holographic 
energy with IR-cutoff on the particle horizon one gets w — 0, which again clearly contradicts the observations. 

Let us convince ourself now that the above discovered problem preserves also i n th e case of Universe dominated by 
holographic dark energy. Present the density of the latter in the following form |l58j | 



PDE = 3c^MpiL-\ (255) 

The coefficient 3c^ (c > 0) is introduced for convenience, and Mpi further stands for reduced Planck mass: Mp^ = 
SttG. Replacing L — >■ Rh in (I255p . one obtains the first Friedmann equation in the form 

RpH = c, (256) 

and it immediately follows that 

1 d / 1 



(257) 



Ha^ da \Ha^ 
One now easily finds that 

H-^ = aa^+-, (258) 

where a is constant. As R^ = cH^^, and pA — "ic^ MpiR^^ , then 

pjjE = 3a2Af|,;a-2(i+l). (259) 

How can we find the parameter w in the state equation if the dependence p = f{a) is given? One possible approach 
is the following. Let us use the conservation equation 



p + iHpil + w) = 0. 
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In the case p = f{a) it takes the form 

f'{a)h + 3Hf{a){l + w) =0, 

and it follows that 

If (a) I dlnfja) 

It is easy to see that for well-known cases of matter p cx a~^, radiation p cx a^"* and cosmological constant p — const 
one obtains the correct values w — 0, 1/3,-1. From the expression (|260p it follows that 

12 1 , , 

- = -3 + 3^>-3- (2^1^ 

Thus one can see that the above described component does not deserves the name of "dark energy" in proper sense, 
as it cannot serve its main purpose — to provide the accelerated expansion of universe. The difficulties originate 
from the fact that, as it follows from (|257p . the derivative obeys j^[H~^ /a) > 0. In order to o btain the accelerated 
expansion of Universe, one should "slow down" the growth of the IR-cutoff scale. It turns out jl23j . that it can be 
done if one replaces the particle horizon by the events one. Recall (see Chapter 2), that the size of spatial region such 
that a signal emitted at time t from any point of it will reach the immobile observer at infinitely far future, equals 

dt' 

Re = a{t) / — -. (262) 
Jt a (f) 



Assuming again that the dark energy dominates and finding solution of the first Friedmann equation in the form 
da c 



(263) 



Ha^ Ha ' 
one finds 

pjjE = ic^MliR-^ = 3a2M|,;a~2(i-i)^ (264) 

orw = — i — -i. We obtained a component which behaves as the dark energy, i.e. it provides the accelerated 

expansion of Universe. If c = 1 then it behaves as the cosmological constant. If c < 1, then w < —1. Such value of w 
responds to phantom model in the traditional approach. 

In the above considered cases we investigated properties of holographic dark energy in two limiting cases: the 
matter dominated and dark energy dominated one. Let us now consider the c ommon case situation, i.e. the Universe 
dynamics at arbitrary relation between the both components densities |124| . For simplicity we restrict to the case 
of flat Universe and set the IR-cutoff scale to the event horizon (35). Introducing the relative density of the 
holographic dark energy VIbe ~ Pde/ Per {Per = •^MpiH'^), we represent (|255p for L ^ Re in the form 

HRe = (265) 

V "-D-E 

Of course, with Ude = ^ and replacement Re Rh the equation (|265p transits into ()256p . Taking derivative with 
respect to time from both sides of (|262p . one gets 

Re = HRe - 1 = , ^ - 1 (266) 
From the definition (j255|) it follows that 



^ = -QcHll,R-^Re = -2H (l ^ j p^E. (267) 
Due to the energy conservation law 

^ {a^PDE) = -ia'pDE- (268) 
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and it follows that 

PDE^-i;-p. Pde- (269) 

6 a ma 

Thus the state equation reads 

Pde I din Pde , 1 



wde 



Pd e 3 d\na 3\ c 



-l = -:7 l + -VfW . (270) 



We used the fact that dlna = Hdt. This result can be obtained without calculation of the pressure pA by the 
relation (|260p . The obtained expression for Wj^^ is a consequence of the holographic dark energy definition in form of 
(|255p . therefore it is independent on other energy components. It follows from the obtained result that w^j^ ~ — 1/3 
in the case of domination of other energy components and w^^^ = — 5 (l + f) in the dark energy dominated case. 
The latter result coincides with the expression (j260p , obtained above for the Universe filled by solely holographic dark 
energy. 

For the first impression the declared task is completed. The holographic dark energy with density (|255l) from the 
one hand provides correspondence between the observed density and the theoretical estimate, and from the other it 
leads to the state equation which is able to generate the accelerated expansion of Universe. However the holographic 
dark energy with IR-cutoff on the event horizon still leaves unsolved problems connected with the causality principle: 
according to the definition of the event horizon the holographic dark energy dynamics depends on future evolution of 
the scale factor. Such dependence is hard to agree with the causality principle. 

In order to seek the way out from that dead end let us address once more the cosmological constant problem, 
which is the huge gap between the theoretical estimates and observed values of the dark energy density. The simplest 
type of the dark energy — the cosmological constant — is attribute to vacuum average of quantized fields and can 
be measured in gravitational experiments. Therefore the cosmological constant problem is the quantum gravity one. 
Although the complete quantum gravity theory is absent, a combination of quantum mechanics with general relativity 
can shed light on the question under discussion. 

From the first days of quantum mechanics the concept of measurements, either real or gedanken, played fundamental 
role for understanding of physical reality. General relativity states that classical physical laws can be verified with 
arbitrary unlimited precision. The above revealed connection between the macroscopic (infrared) and microscopic 
scales dictates the necessity of deeper analysis of the measurement process. The uncertainty relations combined with 
general relativity generates a fundamental space-time scale — the planck length Lp; ~ 10~^'^ cm. 

The existence of a fundamental length scale critically influence the measurement process jl27{ . Assume that the 
fundamental length scale is Lf. As the space-time reference frame must make sense, it must be linked to physical 
bodies. Therefore the fundamental length postulate is equivalent to limitations for possibilities to realize the precise 
reference frames. In terms of light signal experiments this means, for instance, that the time interval required for the 
light signal to propagate from A to b and back, measured by clock in A, is subject to uncontrolled fluctuations. The 
latter should be treated as evidence for the fluctuations of metric, i.e. of the gravitational field. Thus the fundamental 
length postulate is equivalent to gravity field fluctuati ons. 

The existence of quantum fluctuations in the metric [l46l - ll48i Il60| [ directly leads to the following conclusion, related 
to the problem of distance measurements in the Minkowski space: the distance t ^ cannot be measured with precision 
exceeding (l45l | the following 



dt = Pt%^t^/^ (271) 

where /3 is a factor of order of unity. Following jl2l| we can consider the result (|27ip as the relation between the UV 
and IR scales in frames of effective quantum field theory, which correctly describes the entropy features of black holes. 
Indeed, presenting the relation (|252p in terms of length and replacing A — >■ 6t, we now recover (|27ip in holographic 
context. 

The relation (j27ip together with quantum mechanical energy-time uncertainty relation enables us to estimate the 
quantum fluctuations energy density in the Minkowski space-time. According to (I247P we can consider the region of 
volume as composed of cells 6t^ ^ ip^t. Therefore such a cell represents the minimally detectable unit of space-time 
for the scale t. If the age of the chosen region equals t, than it follows from the energy-time uncertainty relations that 



recall that we use the system of units where the light speed equals c = H = 1, so that Lpi = tpi = M 



PI 
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its existence cannot be realized with energy less than ^ t ^. Thus we come to conclusion: if lifetime (age) of some 
spatial region of linear size t equals i, then there exists a minimal cell St^ , with energy that cannot be less than 

Est^^t-\ (272) 

From (|27ip and (|272l) it immediately follows that due to th e energy-time uncertainty principle the energy density of 
quantum fluctuating metric in the Minkowski space equals [l46l ll48L Il60| 



A fact of principal importance is that the dynamical behavior of the metric fluctuations density (|273p coincides with 
the above defined holographic dark energy (|252p . (|255l) . though the expression were derived based on completely 
different physical principles. The holographic dark energy density was obtained on the basis of entropy restrictions 
— the holographic principle, while the metric fluctuations energy density of the Minkowski space is caused only by 
their quantum nature, namely with the uncertainty principle. 

The relation (|273l) allows to introduce an alternative model for holographic dark energy [l60j |. which uses the age 
of Universe T for IR-cutoff scale. In such a model 

P^^^^, (274) 

where n is a free parameter of model, and the number coefficient 3 is introduced for convenience. So defined 
energy density (|274l) with T ~ ■, where i?o is the current value of the Hubble parameter, leads to the ob- 
served value of the dark energy density with the coefficient n value of order of unity. Thus in SCM, where 
ffo 72kmsec-^Mpc-^, Qde ^ 0.73, T ~ 13.7Gyr, one finds that n ~ 1.15. 

Now let us address the crucial question, whether the holographic energy density in form ()274p result in accelerated 
expansion of Universe. For simplicity consider the Universe free of other energy components. In such case the first 
Friedmann equation reads 

The Universe age T in (|274p is linked with scale factor by the relation 

Solution of the equation ()275|) with energy density ([274]) reads 

[n{Hot + a)]". (277) 

The integration constant can be determined from the conditions Uq — I. Evaluating the second order derivative from 
the scale factor, it can be shown that the accelerated expansion of Universe takes place under the condition n > 1. 
We note that the above obtained value n ~ 1.15, corresponding to the observed value of the dark energy density, 
satisfies the above mentioned condition, which can be obtained from the conservation equation for dark energy, easily 
transformed to the following form 

., = -1-^. (278) 

Using ((274t . (|274t and ((2751) we present the expression (I278P in the form 

= -l + |-. (279) 

As was repeatedly mentioned above, the accelerated expansion of Universe requires w < —1/3, which is equivalent to 
the above obtained condition n > 1. Like in the previous model of the holographic dark energy, we transit to more 
general case: the Universe where dark energy coexists with matter with density Such Universe is described by 
the Friedmann equation 

+ (280) 
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Transforming to relative densities f2„i — a/pM^ ^'^"^ ^9 ^ sh^m'^ ^ t"h'^ ' present the Friedmann equation 



([280| in the form 



^"^ = (3_ ^281) 



din a n 
The equation (|28ip can be solved exactly to give 

-\na + co - --^ln(l- yn;)- -^ln(l + v/f^) 
n Sn — 2 6n + 2 



-^\nn,+ ln{^-^^), (282) 



3n * 3n(9n2 - 4) '2 

The integration constant can be determined from the condition flq ~ 0.73 at a = 1. Let us analyze the dynamics of 
relative dark energy density in two limiting cases: the matter-dominated and dark energy-dominated ones. In the 
former case it follows from (|28ip that 

« aa^. (283) 

Fast growth of relative contribution of the dark energy, irrespective of n value, results in the dark energy dominated 
era, when 

Qg^l- C2a-(3"-2)/". (284) 
The state equation for the dark energy can be obtained from the relation (|278p . 

Wq^-l + -^^q. (285) 

In the early Universe during matter-dominated era one gets V,q ^ and Wq —1, i.e. during that period the 
holographic dark energy in the considered model behaves like cosmological constant. In later epoch of dark energy 
domination, when fig — )■ 1, the state equation (j285p naturally transforms to the above obtained relation (j279p . Remark 
that fate of the Universe filled by matter and holographic dark energy with density ()274p is constant accelerated 
expansion with power law time dependence (|277p of the scale factor. Thus the holographic model for dark energy 
with IR-cutoff scale set to the Universe age, allows the following: 

1 . to obtain the observed value of the dark energy density; 

2. provide the accelerated expansion regime on later stages of the Universe evolution; 

3. resolve contradictions with the causality principle. 

However the reader should not hurry with the ultimate conclusions. The first successes of holographic principle 
application from one hand awoke hopes to create on that basis an adequate description of the Universe dynamics, free 
of a number of problems which the traditional approach suffers from. From the other hand, it was those successes 
which provoke we would say unreasonable optimism. We believe that papers entitled somewhat like "Solution of 
dark energy problem" (200l | represent manifestation of specific "holographic extremism" . We would like to remind the 
saying of W.Churchill: "Success means motion from failure to failure, not loosing the enthusiasm". The holographic 
dynamics is one of the youngest directions of theoretical physics. On that path the physicists experienced yet too few 
failures in order to pretend for complete success. 



X. TRANSIENT ACCELERATION 



Unlike the fundamental theories, the physical models only refiect our current understanding of a process or a 
phenomenon, to describe which they were created. The model efficiency significantly depends on its flexibility, i.e. 
its ability to be modernized using the new-coming information. That is why evolution of any actively living model 
includes multiple generalizations, directed both to resolve the conceptual problems and to describe ever growing set 
of observational data. In the case of SCM those generalizations can be divided on two types. The first includes 
replacements of cosmological constant by more complicated dynamical forms of dark energy, taking into account the 
possibility of interaction of the latter with dark matter. Generalizations of the second type are of more radical nature 
and pretend one consequent change of cosmological paradigm. The ultimate (explicit or hidden) aim of them is to 
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reject the dark components due to modifications of Einstein equations, and consequently the Newtonian laws two. The 
generalizations of both types can be well viewed on example of the phenomenon called the "transient acceleration" . 

As we have seen above in frames of SCM, the dependence of the deceleration parameter q on the redshift z has 
the characteristic monotonous trend to the limiting value q{z) — 1 for z — >■ — 1. It physically means that after the 
beginning of dark energy-dominated era (at z ^ 1), the Universe in SCM is condemned to eternal acceleration. 

Below we consider some cosmological models with dynamical forms of dark energy, which lead to transient acceler- 
ation, and we also discuss the information on current Universe expansion rate extracted from observational data. 



A. Theoretical background 



J. Barrow [38[ was one of the first to consider principal possibility of the transient acceleration. He showed that 
many well established scenarios, consistent with the current accelerated expansion of Universe, do not exclude the 
possibility of recurrence to non-relativistic matter-dominated era, and thus to the decelerated expansion regime. 
Therefore the transition to the accelerated expansion does not yet imply the eternal accelerated expansion. 

In order to show that we, following the Barrow work [ssf . consider a homogeneous and isotropic flat Universe, filled 
by non-relativistic matter and scalar field with the potential energy V{(p) and state equation p = wp. We take the 
scalar field potential in the form 

Vi^) = Vp{ip)e-^'^. (286) 

In some versions of low-energy limits of string theory the potential Vp{(p) represents a polynomial. The exponential 
potential with a shallow minimum was first suggested by Albrecht and Skordis (37| . This minimum on the exponential 
potential background was created due to the polynomial factor Vp{(p) of the simplest form 

Vpiifi) ^{ifi- ifio)^ + A. (287) 

In the considered case the potential takes the following form 

Vi^) = e-^^{A+i^-^of). (288) 

In order to relate the above considered potential to that of string theory the constant parameters A and ipo should 
be of order of unity in Planck units. In such quintessence models the accelerated expansion of Universe appears 
naturally on late stages of evolution without any fitting of initial parameters, thus this model is free of the precise 
tuning problem. 

The accelerated expansion starts when the field rolls down the local potential minimum at = <j)Q+{l±\/l -~ X^A) /A, 
which is formed by the quadratic factor in (|288p . where 1 > X^A. While the field stays in the false vacuum state its 
kinetic energy is negligibly small (0 ~ const), and consequent domination is caused by almost constant value of 
the potential energy, which runs the era of the accelerated expansion of Universe which never ends. 

It was discovered in the paper [38^ that it is neither solely possible nor most probable scenario. 

The transient vacuum domination appears in the two following cases. When AX^ < 1, the field (p reaches the local 
minimum with kinetic energy sufficient to overcome the potential barrier and continues roll down the exponential 
part of the potential to the region where ip ^ tpo- The kinetic energy is defined by the scaling regime and thus by 
the parameters of potential rather than by the initial conditions. The transition acceleration appears also in the case 
when the condition AX^ > 1 holds. Since A grows proportional to A~^, the potential looses its local minimum and 
fiattens near the inflection point. It is sufficient to cause the temporal acceleration of the Universe expansion, however 
the field never stops rolling down the potential and once the Universe will be again dominated by matter with the 
dependence a{t) cx t^/^. 

Therefore the Universe quits the regime of eternal accelerated expansion and returns to the decelerated expansion. 
Besides the the well motivated family Albrecht-Skordis potentials the possibility of transient acceleration regime 
appears to be more probable than the eternal accelerated expansion. 



B. Different models with transient acceleration 



In order to show explicitly that the transient acceleration represents a natural feature of different cosmological 
models, we briefiy consider some of them below. 

Barrow considered in his work (ssj a model of Universe where dark energy is present in form of the scalar field with 
potential (|288p . We consider some other examples of cosmological models that provide alteration of the accelerated 
expansion phase by the decelerated one. 
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1. Scalar field, multidimensional cosmology and transient acceleration 

In the paper [66] it was shown that the transition acceleration phase can be reahzed in the exponentail potential as 
well, and besides that a d-dimensional cosmological model was also considered in the work. The action of the latter 
reads 



(289) 



where = SttGn = 1/2 and Vb > 0, A > (the case A < is connected to the case A > by the replacement 
(j) — > —</)). Following the paper [6g| we consider below the FRW-metric for flat Universe (k — 0) 

ds^ = dt^ ~ a^{t)dx'dx\ i = (290) 

In that case the action takes the form 



d 



i'^x (^d ~ l){d ~ 2)a'^-3^2 ^ 0^^-1(1^2 „ v{c^))^, (291) 



and in variables (u, v) 



K\d — 
one obtains 



^ "^(v-u) , a-^-i = 6"+", (292) 



S = f d^-^xdt e^+^ ( ^\ uv - Foe-2"(-'") 

J \K^{d-l) 



Let us transform to new time variable r 



J(^_lMe"(— ). (294) 
dt y 2(d-2) ^ ' 



g^}_ md^Wo I ^d^i^^^ e"+''e"("-"'(uV- 1). 



and therefore 



K V rf — 1 

Using (j293l) . (j294l) and (|292p it can be shown [66], that the general solution for the case a < 1 takes the form: 
2 2{d-2) 4o£.(l + me-2"'^)a^_, 2 4., 2^, , 2 2^, , " , ; 



1 Id- 2 /2ar 



1 . 1 



kVc?— l\w 1 + a 1 — a 

where s — d ~ I, and m is the integration constant. 
The solution (I295P has the asymptotes: 



log(l - me-2"'^) + log(l + 7716"^"'^) ) , (295) 



a^i'i-i , - — ^logt, for t = 
K \ d — 1 

2 



a ^ t(<i-2)A2 ^ = - log t , for t > 1 . (296) 
A 

On early stages of evolution the state equation is extremely rigid p = p, and on later stages it is easy to show that 

d-2 ,2 
P^^'P' "^2.2(d_l)^ -1- 
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In accordance with the with the equation (|296p . the accelerated expansion wih continue forever under condition 

A < J^J^__2 • ^'^^ show that under condition J^J^_^ < A < 2k^ the solution (I295P with m > has a transient 

acceleration stage. First let us find ^ = Using (|295p it is easy to show that a is proportional to positively 

defined quantities m and t. Setting \m\ = 1, which can be easily done by shift along the t axis, one finds a in the 
following form: 

d = -(positive) ^((d - l)^^ - l)Z^ - 2{d - 2) sign(m)aZ + d - 1 - j , (297) 

where Z = cosh(2ti;T). If (d— l)a^ < 1, which corresponds to the case A < ^^"^^ , then we obtain the eternal accelerated 

expansion for arbitrary values of m, because only the first term in ([257]) dominates on the later stages of evolution. 
Remark that such solution represents an attractor. If {d — l)a^ > 1 holds for later times then the solution always 
corresponds to to decelerated expansion of Universe. For {d— l)a^ > 1 and to < the right hand side of the equation 
(j297[) is negatively defined, which corresponds to decelerated expansion of Universe during whole evolution time. At 
last, in the case (d — l)a^ > 1 and to > 0, the solution always provides the transient acceleration phase, which is 
mentioned in[63|. Indeed, the equation (|297p has two roots 



{d~2)a±Vd^{l-a^) 
^± = id-l)a^-l ' 

defined in the interval T-{a) < r < T-i-{a), corresponding to the accelerated expansion. 

The root t± is real (and positive), because on the interval a e (l/Vd — 1, 1) we have Z± > 1. In the limit a —>■ 1 the 
two roots coincide and duration of the accelerated expansion phase tend to zero. In the opposite limit a — > 1/^/d — 1 
one obtains Z_ = d/ {2\/d — 1) and Z+ 00, which corresponds to infinite period of accelerated expansion. For 
higher dimensions of space one obtains Z± — l/a ± (1 — a^)/ {a'^\/d) + 0{l/d), and it means that the transient 
acceleration duration is shorter in the Universe with larger number of spatial dimensions. 

2. Transient acceleration in models with multiple scalar fields 

The existence of an event horizon in the case of a de Sitter phase is an obstacle to the implementation of string 
theory because the S-matrix formulation is no longer possible, hence eternal acceleration, leading to a de Sitter space 
in the asymptotic future is problematic (70| . The transient acceleration models allow to avoid that contradiction, 
because they produce the accelerated expansion today and in recent past, without any contradiction in future. In 
the paper [68| the dynamics of homogeneous and isotropic Universe was considered, where the dark energy is realized 
by two (generally speaking, turn by turn) scalar fields $ and tJj. Let us write down the equations of motion for that 
system 

(p = -SHip^d^V (299) 
^ = -SHijj-d^V, 

and also the first Friedmann equation: 



r2 SttG, , k 
f = ^^[Pb + PQ) 2 

O 0/ 



Here the dot denoted the derivative with respect to time t, the subscript b marks the background components such 
as dark matter (to) or radiation (r), and Q stands for dark energy in form of the two scalar fields. It is convenient 
also to define 

Pb,Q = {lb,Q - l)pb,Q , 

then 7,„ = 1 and 7^ = |, where, for any component i, we have introduced for convenience the quantity 7, = 1 + Wj. 
The energy density and pressure for the quintessence fields in the potential V{p, ip) have the form: 

PQ = + + V{p, i,) PQ = + „ y^^^ 
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As always, the possibility of the transient acceleration crucially depends on the potential form. Below we consider 
several possibilities to obtain the transient acceleration. The paper [68j considers a number of examples, including 
the case when the two fields interact through the coupling potential V{^p,^lJ), as well as the case when one of the 
fields is free. Here we consider only the models with potentials depending on both fields. Direct generalization of the 
potential from the paper (I288P for the case of minimal coupling between the two scalar fields enables us to obtain the 
transient acceleration. 

Consider a potential of the form 

V{^, ^) = M^e-^'^{Po + - ^c)' + gm- (300) 

The above introduced additional scalar field ip will control absence or presence o the potential minimum for ip. The 
main idea is that, from one hand, a potential minimum initially exists for the scalar ip, which provides accelerated 
expansion of Universe, and from the other hand evolution of the field ip results in the extinction of the minimum and 



the Universe returns to the decelerated expansion regime. The model of [37| is recovered with / = 1 and g = 0. For 



the potential (|300[) the position of the minimum is determined by the expression 




The function g {g > 0) is responsible for the scalar field mass and it usually takes the form g oc ip"^, but as it does not 
affect the system dynamics, we set for simplicity g = 0. 
The minimum (|30ip disappears under the condition 

fiij) < A^Po = fiyJc). (302) 

Note that the potential pOOp can be presented in the form 

V{^,^) = M^/\^e-^^{f{iP,) + fm\^ - \^,f). 

We assume that / is positive, continuous and monotonous function for -0 > and/or -0 < 0. From the condition (|340|) 
it is easy to see that if / decreases (increases), then for ip < tpi {ipi— is the initial value of the scalar field), {ip < ipc) 
the accelerated expansion takes place, because the minimum poip exists. One is free to use different functions / to 
obtain the transient acceleration, for instance we present the case p, a > 0, / = 1 + aijjP, f = ta,iih{aipP), or for 
arbitrary p and a, / = exp{a'ipP) , f = coshlatpP). During the evolution of Universe the scalar field (p rolls down the 
potential and dominates in the exponential part so that Mpim^ ^ Mpiirn^ V ^ tf^ ^ Mp^H"^ with Q,q = 
wq = 1/3 during the radiation dominated era, and with 51q = 3/A^, wq — for the matter-dominated one, while ip 
approach the minimum. When the field p is in the minimum ()301|) and oscillates near its minimum value, the universe 
expands with acceleration with V ^ pP' and therefore wq ~ — 1 until ip obeys < 'i/'c, ("0 ^ V'c), if the function / 
grows (decays). In that moment the minimum (|30ip disappears and the field ip starts free roll acquiring large values 
(speeding up) and the the Universe expansion becomes decelerated one q > Q. 

When -0 is initially larger (smaller) than 0c, provided / is increasing (decreasing) passes through t/'c because 
d^V — M^e^^'^{(p ~ (pcY^ is positive (negative), acceleration occurs which is always transient; if on the contrary 

"04 ^ V'c {''Pi > '0c) quintessence domination is not possible. The critical value ipc controls the presence or not of the 
minimum for 0. 

For a particular example let us consider the simplest case 

f{iP) = iP^ , (303) 

then the minimum for p disappears under condition —ipc < ip < ipc, where ipc = \^/Pq. In analogy with the model 
()288p . one can obtain A > 9, which agrees with the limitations imposed by the cosmological observations, provided 
the initial value p)i is fixed and the critical value ipc must be finely tuned in order to make the quintessence dominate 
in the present time. 

The diagram |9] shows the state equation parameter wq^ ~ —0.491 and Weff,o = —0.874. When the field (p gets into 
the minimum, then wq ~ — 1. The deceleration parameter q is also shown on the plot, and it points out on the fact 
that in the present time the Universe already passed the transient acceleration stage and currently it expands with 
deceleration qo ~ 0.013. The accelerated expansion {q < 0) starts at z ~ 0.658 and stops at z ~ 0.0035, when the 
Universe age equals tend/to — 0.996, while in the present time one has Hoto ~ 0.912. 
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FIG. 9: The diagram plots the dependence of the state equation parameter lo^ (solid line) and deceleration parameter q (dashed 
line) as functions of the redshift z for the potential pOO[l with the fixed parameter values A = 10, </>c = 23.8 and Po taking the 
following values values (upward from below) 0.160,0.162,0.163,0.164,0.166,0.168,0.170 {4>i = and ipi = 5). In the models 
with Po ~ 0.164, 0.166, 0.168, 0.170 the accelerated expansion is replaced by the decelerated one in the modern epoch. 



Decaying dark energy as a scalar field 



As was already mentioned above, the cosmological models where the dark energy decays have many attractive 
features, one of which being the presence of transient acceleration — alteration of the accelerated Universe expansion 
period by the decelerated one, which takes place when the dark energy density becomes sufficiently low. Dynamics 
of such Universe clearly differs from its evolution predicted by SCM. The considered model with the decaying dark 
energy represents a pre-image of the commonly accepted inflation model, where the field, generated the inflationary 
expansion of Universe, experiences the decay. 

Consider a scalar field model with the potential that takes both positive and negative values 



107|: 



(304) 



The assumption that the scalar field potential can take negative values is very curious, but one should remember that 
except the motivation to explain the observed transient acceleration this model lacks any other observational support. 
Nevertheless such effective potentials often appear in the super-gravity and M-theory. It is be shown below how to 
obtain the transient acceleration in frames of that model. Remark also that evolution of Universe in cosmological 
models with negative potentials sharply differs from that of SCM. So, for example, in the case when the Friedmann 
equation is dominated by potential energy p04|) . taking negative values 



1 



1 



(305) 



even the spatially flat Universe can collapse, which is principally impossible in SCM, but it can happen only on 
later stages of Universe evolution, far before the transient acceleration ends. Up to the moment when the Universe 
starts collapsing, many stages of transient acceleration manage to occur (see Fig llOp . which become more often when 
approaching the collapse moment. 



57 



2,0- 



1,5- 



1,0- 



q 



0,5-- 



0,0 



deceleration 



acceleration 



-0,5- 



-1,0- 



2 3 

a 



FIG. 10: Dependence of the deceleration parameter on the scale factor for the Universe filled by matter and dark energy in 
form of scalar field in the potential V{Lp) = Voces (ipm/y'V^) , with m = 0.74; Vq = 150, 92(0) = 0.23, 93'(0) = 0. The present 
time corresponds to a = 1. The potential parameters are chosen so that the deceleration parameter value appears to be the 
same as in SCM, q{l) ^ —0.6. 



4- Transient Acceleration In the Universe with the Interacting Components 



Consider the spatially flat Universe composed from three components [101|: dark energy, dark matter and baryons. 
The first friedmann equation then takes the form: 



?>H^ = p 



DE ' Hrn 



Pb, 



(306) 



where as usual is the dark energy density, p„ ~ dark matter energy density, pi, - the density of baryons, H — ^ 
- Hubble parameter. The state equation for the dark energy takes the form = wp^,^ . 
The conservation equations in that case take the form: 



PnE +3i?(i + = -Q; 

Prn + iHpra = Q, 

where Q characterizes the interaction. The conservation equation for the baryons reads: 

, 3 



(307) 



(308) 



Total density equals p = pm + Pb + Poe ■ Without lack of generality, we assume that the dark matter energy density 
satisfies 



Pm=Pmo(^) /(a): 



(309) 



where pmo and oq are constants and /(a) is arbitrary time dependent function. From p07p and p09|) one obtains 

(310) 



Q = Pmj = PtuO y — j J- 



Let /(a) takes the form 
f{a) = l+g{a). 



(311) 
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Interaction free case corresponds to /(a) = 1, thus the function g(a) responds for the interaction. Taking into account 
that 

/ = .^=|«, (312) 
one obtains 

Q = /5m0^a(-^) • (313) 
For pm it means that 

Pm= Pmoil + 9)(—)\ (314) 



where Pmo = Pmio-o) in presence of interaction, and analogously pmo in uncoupled case. The two initial values of dark 
matter density are linked by the following relation 



PmO = PmO (1 + go) , (315) 

where go = g{ao). As can be seen from (I310p . in the case when Q is positive the dark energy decays into dark matter 
^ > and vice versa ^ 



^ > and vice versa < otherwise. It follows from the equations p07p and pi3p 



PnE+^H{l + w)p^^ =~p^o^a(—) . (316) 



da \ a 

Assuming that w — const the solution of equation pi6p takes the form 

'ao\ 3(1+1") 



(Pdeo + Pmogo) (^) 



PoE 

3 ' " ' .a (317) 

'P,nO (^) .9 + Swpraoala~^^'+^^ I daga^^-' 



Rewrite the second Friedmann equation in the terms of g{a) 
! = -i {-Prno (1 + 9) i^f + pm {jf + (1 + 3Z.) X 

X \poEo+PM){jf'^'^'"^ -p™o(^)% + 3z«p™oaga-3(i+»)/;_daga3-i]}. 

In order to solve the equation piSp one have to explicitly define the function g{a). For the lack of knowledge about 
nature of the dark energy, as well as the dark matter, the function 17(0) cannot be derived from the first principles. 
Therefore for the model under consideration we introduce the interaction such that the resulting dynamics in the 
model corresponds to the observed data. 

Consider interaction with the function 17(0) in the form g (a) — a" exp (— a^/cr^) , where n is natural number, and 
a is positive real number. The transient acceleration implies that the dark energy density starts to diminish, i.e. 
experiences decay ^ > 0. This condition requires that n and a satisfy the inequality ncr^ > 2. The dependencies of 
relative densities on the scale factor for n — 7 and a = 1.5 are shown on Fig. 1111 The considered model which allows 
to produce the transient acceleration with certain choice of the interaction parameters, however for large values of 
scale factor as well for small ones it is indistinguishable from SCM. 



5. Decaying cosmological constant and transient acceleration phase 



As a simple example of the transient acceleration we consider a model [lOlj with decaying cosmological constant: 



Pdm + i-Pdra = "PA , (319) 

a 

where /9(j„j and px are energy densities of dark matter and cosmological constant A respectively. On early stages of 
Universe expansion, when is sufficiently small, such decay has almost no influence on the cosmic evolution. On 
later stages, the more the dark energy contribution grows, the more its decay affects the standard dependence for the 
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FIG. 11: The dependencies of relative densities on the scale factor for n = 7 and g — 1.5 (on the left). On the right: the 
dependencies of deceleration parameter on the scale factor in the model with interacting dark energy and dark matter g(a) (solid 
line) with n = 7 and a = 1.5 in comparison to SCM (dashed line). 



dark matter energy density pdm c>c a on the scale factor a. Let us assume that such deviation can be described by 
a scale factor function e(a), then 

Pdm = Pdm,oa^'+'^°^ , (320) 

where in the present epoch gq — 1. Other material fields (radiation and baryons) evolve independently and conserve. 
Therefore the dark energy density takes the form: 

/■^ e(a) + &e'ln(&) ^ 

PA = Pdmfl J a^-iiia) + ^^^^^ 

where prime denotes the derivative with respect to the scale factor, and X - is the integration constant. Neglecting 
the contribution of radiation, the first Friedmann equation takes the form 

H^Ho [f)6,oa-' + ^drnMa) + ^x,q] . (322) 
The function <p(a) reads 

^(a) = + '-^^^^^d,, (323) 

where flx,o denotes relative contribution of the constant X into the total relative density. In order to proceed further, 
some assumptions on the explicit form of the function e{a) should be done. In the present review we follow the original 
paper [101] on A-dark matter interaction and consider the simplest case 

e(a) = eoa^; 

= eo(l + ^)-«, (324) 
where eo and £, can take both positive and negative values. From the above cited expression (|32ip it follows that 

PA = Pmoeo C ^\tY^fd} '^~^ + ^- (325) 

Remark that the case eo = corresponds to SCM, i.e. X = paq. 

Using the above cited formulae it is easy to obtain the dependence of the relative densities VL]y{a), fldmio-) and 
fiA(a): 

- A + a-3 + B-V(«) ^ ^'"'^^^ 
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a-3+£(a) 

= n^R -3^ ( V (326b) 

D + .^(g) - a-3+^(°) 
D + Ba 3 + (^(^aj 

where A = nx,o/^b,o, B = rib^o/^^dm.o and D = f^x,o/f^dm,o- 

In such simple model, appropriate choice of the parameters eo and ^ enables to obtain practically any kind of 
Universe dynamics. In the context of the present review it is especially interesting to consider the case when eo > 
and ^ takes large positive values > 0.8). The figure [T^ shows the dependence of the scale factor for the case ^ = 1.0 
and eo = 0.1. Note that with such parameters in the present time a 1 the Universe expands with acceleration, but 
unlike the case of SCM, the dark energy domination will not last forever and at a ^ 1 the Universe enters into new 
era of non-relativistic matter domination with a — ^ cx). Such type of dynamical behavior is impossible for most models 
with A{t) or those with interacting quintessence, widely discussed in literature, but it is an attribute of the so-called 
thawing '4^ and hybrid [4TI| potentials, which follow from string or M-theory (see also [43j)^. 

In order to give clearer idea of transient acceleration phenomenon we find the explicit form of the deceleration 
parameter q = — aa/a^, in the considered model 

. ^ 3 r!fc,oa-^ + »d„,oa^(°)-3 

Its dependence on log(a) for selected values of ^ and eo is shown on FigH^l Note that large values of the parameter 
^ corresponds to the Universe dominated by matter in the past {q{a) — >■ 1/2 for a 3> 1), afterwards (at Qacc < 1) 
a long-during era of accelerated expansion takes place, which, unlike SCM, does not last forever, ending at some 

a-dec > 1- 



C. Transient acceleration: holographic limitations 

The novel concept of the holographic Universe considered in the previous section turns out to impose strict lim- 
itations to Universe dynamics. In order to see that we consider the matter entropy contained within the observed 



Authors of the paper [S^l present arguments for the fact that forever expanding Universe, which is common feature of most quintessence 
models (including the standard ACDM- model), contradicts the predictions of string/M-theories, because those models contain cosmo- 
logical horizon, which make impossible to apply the usual procedure of S-matrix formalism, describing the interaction of the particles. 
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horizon in the case when the matter is composed of non-relativistic substance. Each particle of the substance carries 
one unit of information, and therefore Sm = KqN, where N — (47r/3)rJ^n, and n is the particle concentration, for 
which n = noa~^. Thus one obtains for entropy the following 

4:7r 

S^^ks — f\ no cx a^/^ ^ (328) 

and it follows that S'J^ oc > 0. It is interesting to note that as the Universe expands the particle number contained 

inside the horizon increases too, which provides the entropy growth and therefore S'J^ + S'^ > 0. It is easy to see that 
the Universe tends to maximum entropy provided its horizon decreases with time. In the simplest case the horizon 
equals to Hubble radius and the above condition is satisfied if the Universe expands with acceleration. In such case 
the deceleration parameter is negative {q < 0), thus 

^{Rh) = c{l + q) < c, 
at 

and Hubble sphere has velocity which is less then light speed by the quantity cq, and thus it is left behind those 
galaxies. Therefore the particles initially placed inside the Hubble sphere, gradually leave it. 

From the latter fact that a conclusion follows that in order to obey the second law of thermodynamics, the universe 
must be dominated by dark energy with state equation parameter satisfying the condition —1 < w < —2/3. Another 
alternative is to modify the gravitation theory in order to obtain the accelerated dynamics without the dark energy. 

Likewise, since every isolated system is expected to evolve in a such a way that it tends to thermodynamic equilib- 
rium the curvature of the A{a) function should never spontaneously change from negative to positive values. In order 
to specify the set of possible cosmological values, or equivalently to impose certain limitations on the values taken by 
the state equation parameter w, it is necessary to use the laws which are not immediately connected to cosmology. In 
the context of holographic cosmology it is reasonable to use the generalized law of the black hole thermodynamics — 
the generalized square law (GSL). According to the latter the total entropy of the system S must obey the condition 
S" < 0. As will be shown below, (see jl04{ ) application of this criterion on later stages of Universe evolution helps 
considerable specify the range of values allowable for the state equation parameter w. 

For further convenience we express the horizon area A and its derivatives through the deceleration parameter 
q = —d/{aH'^), which can be presented in the above mentioned form: 

.= -(l + ^). (329) 
For the case of spatially flat Universe A — AttH^^, therefore A' — 8tt{1 + q)/ H"^ and 

A!' ^2A (1 + q) (1 + 2a) + ^ . (330) 



(1 + q)[l + 2q) + ^ 
a 



It follows that for the case q > —1/2 the second order derivative of the horizon area. A", cannot change from positive 
to negative values while q' < 0. It excludes the dependence for q shown on FigH^l Such dependence of the deceleration 
parameter appears in the cosmological models where the current accelerated expansion is transient, and after a short 
dark matter d ominated period the dark matter starts dominating again, which leads to the decelerated Universe 
expansion era. [lOSl Il06j | . 



D. Transient acceleration in models with holographic dark energy 



Current literature usually considers the models where the required dynamics of Universe is provided by one or 
another, and always only one, type of dark energy. As was multiply mentioned above, in order to explain the observed 
dynamics of Universe, the action for gravitational field is commonly complemented, besides the conventional matter 
fields (both matter and baryon), by either the cosmological constant, which plays role of physical vacuum in SCM, or 
more complicated dynamical objects — scalar fields, /^-essence and so on. In the context of holographic cosmology, 
the latter term is usually neglected, restricting to contribution of the boundary terms. Nevertheless such restriction 
has no theoretical motivation. 

In the present subsection we consider the cosmological model which contains both volume and surface terms |l78l 

Il79l |. The role of former is played by homogeneous scalar field in exponential potential, which interacts with dark 
matter. The boundary term responds to holographic dark energy in form of (|28ip . Such model is shown to produce 
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a 

FIG. 13: The deceleration parameter dependence in some models with transient acceleration. 



the transient acceleration phase — the evolution stage when the accelerated expansion of Universe changes to the 
decelerated one, and afterwards the Universe comes to the stage of eternal accelerated expansion, thus removing the 
limitations imposed in the previous subsection. 

To describe the dynamical properties of the Universe it is convenient to to transform to dimensionless variables as 
the following: 



y = T7-7rV^^. ^ = T7-TT\hr^ ^^T7-TT\hr- (331) 



VGMpiH' ^ AIpiH\ 3 ' MpiH\ 3 ' MpiH \ 3 

The evolution of scalar field is described by the Klein-Gordon equation, which in the case of interaction between the 
scalar field and matter takes the following form: 

!f + 3Hip+ — ^~-^. (332) 

dip ip 

In the present section we consider the case when the interaction parameter Q is a linear combination of energy density 
for scalar field and dark energy 

Q^3H{ap^ + (3p^), (333) 

where a, /? are constant parameter. For given model, regardless the explicit form of the scalar field potential V{p)^ 
the system of dynamical equations takes the following form 

3x , , „ / 3 , 2 



where 



and 



y5(x,z,M) -3a;+ y-Ay -7, 
y' = Y5(a;,2,u) - y^Axj/, (334) 

w = -:r9{x,z,u) , 

2 n 

g{x, z, u) = 2.t2 + z^ + A^3^ ^ ^ _L^Mpi. (335) 
6n V dp 

7 = <-'+y')^P^\ (336) 
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As was mentioned above, here we consider the simplest case of exponential potential 

' 2 /i(y3 \ 



V = Vq exp 



3M 



PI 



(337) 



where fi is constant. Taking into account the expression p36p . the system of equations (I334p reads 

a{x'^ + y2) + 



g{x,z,u) 



3x 
Y 



—g{x,z,u) 



l{x,z,u) + 



3z 
~2 



3it , N 

—g(x,z,u) . 

2 n 



[ixy. 



3x — fiy 



— z, 
2 ' 



(338) 



In the considered model the deceleration parameter takes the form 

2 



2a:^ + H 

3n 



(339) 



Note that the cosmological parameters have no explicit dependence on the interaction form, and only determine the 
evolution dynamical variables. This fact essentially complicates analysis of such system. Let us make few remarks on 
the system of equations ([335]). First consider the case y = 0, which corresponds to free scalar field. It easy to obtain 
some restricting relations on the interaction parameters, which follow from the requirement for energy density to be 
real 



2\l- > l + a + (3, 
a 



(340) 



which in turn requires > /3 > a, |a| + |/3| < 1. This unstable critical point corresponds to Universe filled by dark 
matter. The case of multiple such points is also possible but irrelevant. To conclude we note that any of such critical 
points is easily shown to exist also in the interval xq < 0. For the case z 7^ with the restrictions imposed by p40|) . 
one obtains 



1/2 



where a = - (1 + a + /3)^ 



(341) 



1/4 



1. Case Q = 



In the present section we consider in more details the interaction-free case with scalar field and dark matter. The 
critical points of the system (I338P in the case (a = /? = 0) are given in the table ID The phase space generated by the 
system of equations (|338p contains six physically relevant critical points, the last of which is an attractor. The first 
critical point (1, 0, 0, 0) is unstable and corresponds to Universe dominated by scalar field with extremely rigid state 
equation (w^p — 1), the second critical point is also unstable and responds to the evolution period when the scalar 
field is dynamically equivalent to cosmological constant. 

The next third point (0, 0, 1, 0) is of no interest because it corresponds to Universe composed solely of dark matter 
which is also unstable. The fourth critical point (0, 0, 0, 1) corresponds to Universe solely composed of holographic 
dark energy in form of (|28ip . and was considered above in details. The most physical interest is presented by the last 
sixth critical point which is an attractor. It corresponds to the Universe filled by scalar field and holographic dark 
energy. This critical point is completely determined by the scalar field potential parameter /x and magnitude of n : 



V V ^ (342) 

z. ^ 0, = ^ -i + ^i + l^iil . 
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TABLE I: Critical points for the autonomous system of equations (|338p 



(Xc,yc,Zc,Uc) 


Stability 


q w^p 


Wtot 


coordinates 


character 






(1,0,0,0) 


unstable 


2 1 


1 


(0,1,0,0) 


unstable 


-1 -1 


-1 


(0,0,1,0) 


unstable 







(0,0,0,1) 


stable 




-1 + f 

3n 


( 3_ _3_ /i 3 


- , 0) unstable 







(a;.,y.,0,?i.) 


attractor 


< w^-t 





1.00 



0.25- 



0.00 




_2 Q 2 '^■^ ''■^ ^■'^ 

Ina z 

FIG. 14: Behavior of fl^p (dot line), fi, (dash line) and (solid line) as a function of Af = Ina for n — 3, a — 0.005 and 
/X = — 5 (left side). Evolution of deceleration parameter for this model (right side). 

The obtained property typical for the so-called tracking solutions ,20 li] . Remark also there is also so-called 

background interaction between the scalar field and dark energy, caused by the fact that the scalar field dynamics is 
affected by the holographic dark energy which has negative pressure and influences on the Universe expansion rate by 
means of Hubble parameter, which enters into the Klein-Gordon equation for the scalar field. In the attractor point 
the dark matter density turns to zero. To fit that model with the observations one has to strictly define the initial 
conditions such that the accelerated expansion of Universe starts before then the considered evolution phase begins. 

2. Review of the case Q — SHap^p 

In the above case, the phenomenon of transient acceleration that occurs in such a Universe does not match the 
observations. In order to fix that problem we consider a model where scalar field interacts with dark matter. In the 
present section we consider the case with the interaction parameter of the form (|333p with /3 = 0. The Figure [HI 
shows the dependencies fig, rim and Vl^p for the case when a = 0.005, fi = —5 and n = 3. From the explicit form of 
the equations and interaction character it is easy to see that neither type nor position of the above obtained critical 
points change when the interaction disappears. The latter affects only behavior of the dynamical variables, which 
corresponds to different trajectories in the phase space between the critical points. This corresponds to the fact that 
the interaction parameters enter only in the Hubble parameter derivatives of second order and higher with respect to 
time. 

With such values of interaction parameters the transient acceleration begins almost in the current era. Like in the 
conventional cosmological models with dark energy in form of scalar field, the latter starts to dominate and causes 
the accelerated universe expansion phase. Along with the Universe expansion the contribution of flq increases, which 
leads to the fact that the background (space) starts increasing faster than the field and becomes asymptotically free. 
Such field is known to have the so-called super-strict state equation and forces the Universe to decrease its expansion 
rate. Soon however, when the contribution flq increases enough that the scalar field cannot any more prevent the 
expansion. Universe starts to speed up again. 
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FIG. 15: The deceleration parameter dependence q{z) reconstructed from independent observational data, including the 
brightness curves for SNe la, cosmic microwave background temperature anisotropy and baryon a cous tic oscillations (BAO), 
SN and parametrization (|344|l . The red solid line shows the best fit on the confidence level la CL [l02l ]. 



E. Observational evidence 

Starobinsky [l02| with co-authors, based on independent observational data, including the brightness curves for SNe 
la, cosmic microwave background temperature anisotropy and baryon acoustic oscillations (BAO), were able to show 
(see FigllSp. that the acceleration of Universe expansion reached its maximum value and now decreases. In terms 
of the deceleration parameter it means that the latter reached its minimum value and started to increase. Thus the 
main result of the analysis is the following: SCM is not unique though the simplest explanation of the observational 
data, and the accelerated expansion of Universe presently dominated by dark energy is just a transient phenomenon. 
The term "transient accele ratio n" is increasingly often used with respect to the Universe dynamics. 

Remark, that the paper jl02l | also showed that with the use of CPL parametrization 

wiz)^wo + ^. (343) 

for the state equation parameter, it is impossible avoid contradictions in attempt to combine the data obtained from 
observations of near supernovae of type SNla with those of the relict radiation anisotropy. One of the way out of that 
contradiction is to reject the above mentioned parametrization and invent a new one. Starobinsky group suggested a 
novel parametrization which can combine those data sets: 

w{z) = _l + tanh[(.-zOA]^ ^3^^^ 

Within such approximation w = — 1 on early times of Universe evolution and increases up to its maximum value 
w ^ for small values of z. 

Figure [15] shows the dependence of the deceleration parameter q reconstructed by the parametrization p44p . 

In the year 2 010 i n frames of Supernova Cosmology Project (SCP), the most recent data set on supernovae outbursts 
was published jl03j |. It contains 557 events and is largest up to date. Besides that the data set concerning supernovae 
with moderate red shift (z < 0.3) was remarkably extended. Today there are several works (ssl. [soj. making analysis 
of the above mentioned data in order to test the transient acceleration hypothesis. All the authors share common 
opinion that the ultimate answer will be given only by remade more accurate measurements. Moreover, obtain a 
convincing result one will probably have to remake all the data treatment procedure. Thus it was shown in jSSl . [soj 
that there is a contradiction between the data obtained from (SNe la+BAO) for small redshifts and those from 
(CMB) for high ones. The contradiction is in the fact that the analysis of two data sets separately one obtains 
opposite results. For example, using only (SNe la+BAO) data, one obtains evidence in support of the fact that the 
Universe expansion rate reached the maximum value at (z ^ 0.3) and presently starts to decrease. At the same 
time, if the data are complemented by the CMB data set, the results of analysis drastically change and no deviations 
from ACDM can be revealed. Therefore the result of reconstruction the evolutionary dependence of dark energy and 
answer to the question whether our Universe will expand with deceleration or accelerated expansion will continue 
forever (like in SCM), strongly depends on data obtained from SNe la observation, its quality, the particular method 
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1.0^ 




0.0 0.5 1.0 1.5 2.0 -2.0 -1.5 -1.0 -0.5 0.0 

FIG. 16: On the left: the deceleration parameter dependence q{z) reconstructed from results of analysis of Union24-BAO data 
sets with confidence level 2a. The grey region and the one between the two dashed fines corresponds to presence and absence 
of systematic errors in the SNe la observations. On the right: the confident regions for 68.3% and 95% levels are shown for 
wo and wi in the parametrization CPL with w — wo + wiz/{l + z). Dashed, solid, thin solid lines corresponds to the results 
of Union2S, Union2S+BAO and Union2S+BAO+CMB respectively. The point wo = —1, wi = corresponds to spatially flat 
ACDM model(SCM). 



of reconstruction for the cosmological parameters such as q{z), w{z) and ^^de), as well as on explicit form of the state 
equation parametrization. For detailed answer to that question we have to wait more precise observational data and 
seek less model-dependent ways of their analysis. 
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